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I.  INTRODUCTION 


Studies  of  the  interior  ballistics  of  guns  and  rockets  and  the  dynamics  of 
weapon  systems  require  the  recording  and  analysis  of  large  amounts  of  time- 
correlated  information.  The  type  of  information  associated  with  ballistic 
phenomena  usually  has  a  high  frequency  content  and  a  short  duration.  Wide¬ 
band,  frequency-modulated  (FM)  magnetic  tape  recorders  are  used  to  record  this 
type  of  data.  These  recorders  have  the  capability  of  recording  simultaneously 
fourteen  channels  of  80  kilohertz  (kHz)  data  at  their  maximum  recording 
speed.  A  problem  arises  with  this  type  of  recording  when  an  automated  data 
reduction  system  is  used.  The  problem  is  that  automated  data  reduction 
methods  invariably  use  a  digital  computer  which  requires  that  the  information 
be  in  binary  format.  The  obvious  solution  is  to  sample  the  original  FM  tape 
at  an  appropriate  rate  and  to  convert  the  information  from  analog  form  to 
digital  form  (A/D)  which  can  then  be  recorded  on  a  digital  tape. 

Since  the  digitizing  equipment  used  in  the  Interior  Ballistics  Division 
( IBD)  of  the  Ballistic  Research  Laboratory  (BRL)  either  was  not  computer- 
controlled  or  has  severe  main  memory  limitations,  the  IBD  recently  purchased  a 
Hewlett-Packard  (HP)  minicomputer,  Model  1000-F,  with  which  to  digitize  firing 
data,  off-line,  from  analog  tapes.  The  biggest  advantage  to  this  system  is 
that  the  available  main  memory  does  not  limit  the  amount  of  data  which  can  be 
digitized  in  one  pass;  this  is  accomplished  using  a  double  buffering  technique 
which  effectively  allows  the  transfer  of  data  from  tape  directly  to  disc.  The 
capacity  of  the  disc  now  becomes  the  limiting  factor  of  the  amount  of  data 
digitized.  Reports  on  the  actual  A/D  equipment  and  the  associated  software 
will  be  published  soon. 

The  data,  once  digitized  and  stored  on  the  disc,  can  then  be  transported 
to  the  mainframe.  At  BRL,  this  is  a  CDC  CYBER  System:  1)  mainframe  A  (MFA) , 
a  CYBER  170/173,  and  2)  mainframe  Z  (MFZ) ,  a  CYBER  70/76.  A  CalCorap  plotter, 
Model  1055,  is  a  peripheral  to  that  system. 

A  procedure  is  described  in  this  report  for  the  step-by-step  processing 
of  these  data  until  they  are  converted  into  engineering  units  and  then  stored 
on  a  permanent  file,  integrated,  tabulated,  and/or  plotted.  The  computer 
programs  involved  are  given  as  well  as  sample  inputs  and  outputs. 

II.  CONVERSION  OF  DATA  ON  ANALOG  TAPE  TO  DIGITAL  TAPE 

An  oscillograph  record  of  a  single  channel  of  interior  ballistic  data  is 
shown  in  Figure  1.  It  is  an  analog  pressure-time  record,  with  calibrations, 
reproduced  from  one  channel  of  a  14-channel  magnetic  tape.  To  convert  this 
record  to  digital  form,  an  A/D  converter  samples  the  analog  signal  at  a 
selected  time  and  at  an  appropriate  rate.  It  is  important  that  the  sampling 
times  and  rates  preserve  the  original  pressure-time  relationship;  i.e.,  the 
sampling  times  must  be  directly  related  to  the  real  time  of  the  event. 

The  sampling  rate  for  the  calibration  steps  is  somewhat  arbitrary.  The 
only  basic  requirement  is  to  get  enough  samples  to  make  a  reasonable 
evaluation  of  the  level  by  averaging.  The  sampling  rates  for  the  record 
portion  must  next  be  determined.  In  this  case,  the  sampling  rate  is  more 
critical.  If  the  rate  is  too  low,  there  will  be  insufficient  information  to 
reconstruct  the  continuous  record.  For  example,  in  Figure  1,  if  a  5  kHz 
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Figure  1.  Analog  Pressure-Time  Record 

sampling  rate  is  used  for  the  pressure-time  portion,  it  would  be  impossible  to 
make  an  accurate  peak  pressure  measurement  from  the  few  samples  obtained. 
Conversely,  if  the  sampling  rate  is  higher  than  that  necessary  to  reconstruct 
the  signal,  expensive  computer  time  will  be  wasted  processing  unnecessary 
data.  It  is  obvious  that  a  million  samples  would  not  necessarily  produce  a 
better  defined  record,  indicating  that  there  is  a  direct  relation  between  the 
rate  that  the  data  vary  and  the  number  of  the  samples  needed  to  accurately 
reconstruct  it.  The  rate  that  a  signal  varies  is  directly  related  to  its 
maximum  frequency  content  or  bandwidth.  The  commonly  used  sampling  rate  is 
five  times  the  maximum  data  frequency  content  or  bandwidth  of  the  signal. 

The  bandwidth  of  an  FM  tape  recorder  is  directly  proportional  to  the  tape 
speed.  The  actual  bandwidth  versus  speed  varies  with  the  recorder;  for 
instance,  using  FM  Wide  Band  I  at  120  inches  per  second  (ips)  the  FM  bandwidth 
is  80  kHz,  at  60  ips  the  bandwidth  is  40  kHz,  etc.  A  proper  sampling  rate  can 
therefore  be  selected  based  on  the  tape  speed.  For  example,  to  record 
information  with  a  frequency  content  of  40  kHz,  the  tape  speed  would  be  set  at 
60  ips  which  would  indicate  a  200k  samples  per  second  rate  that  is  five  times 
the  bandwidth. 

For  data  channels  with  calibration  steps,  one  has  a  choice  of  two 
digitizing  procedures.  Normally  the  calibration  duration  is  approximately  300 
milliseconds,  while  the  event  lasts  only  ten  or  twenty  milliseconds.  The 
first  choice  is  to  make  two  passes  over  the  data:  one  at  a  low  digitizing 
rate  for  calibration  steps  and  one  at  a  high  rate  for  the  event.  Often  this 
method  leads  to  bookkeeping  uncertainties  in  matching  calibration  steps  to  the 
event.  The  second,  and  more  preferred,  choice  is  to  digitize  both  calibration 
steps  and  the  event  at  the  higher  digitizing  rate  and  to  eliminate  the 
uncertainty,  even  though  many  more  data  are  generated.  The  excess  data  can 
then  be  eliminated  at  a  later  time  as  discussed  in  Step  3. 

The  HP  digitizer  has  the  capability  to  digitize  up  to  16  channels  of  data 
in  one  pass  of  the  analog  tape.  Consequently,  the  data  that  are  stored  on  the 
disc  are  interwoven;  that  is,  they  are  stored  Xj,  yp...,  z,,  X2>  Yo****) 
Z2>...  These  data  need  to  be  sorted  and  this  is  discussed  in  Step  2. 
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Next,  the  data  on  the  disc  file  of  the  HP1000  must  be  transported  to  a 
disc  file  on  the  BRL  mainframe  where  the  data  analysis  programs  reside.  The 
data  can  be  transported  in  two  ways:  1)  They  can  be  transferred  to  a  9-track 
digital  tape  which  is  then  hand-carried  to  the  CYBER  tape  library.  2)  They 
can  be  transmitted  interactively  from  the  disc  through  an  HP9845  desktop 
computer/terminal  to  the  CYBER  MFA  via  a  modem  link.  This  process  is 
discussed  in  Step  1. 

Step  1.  Creating  a  file  on  CYBER  MFA  or  MFZ 

Data  can  be  transferred  between  the  HP1000  and  the  CYBER  using  the  modem  link 
on  the  HP9845  or  using  magnetic  tape.  Each  method  is  discussed  below. 

A.  Modem  link 

1.  We  have  chosen  not  to  link  the  HP1000  directly  to  the  CYBER. 

Instead,  the  connection  to  the  CYBER  is  through  an  HP9845  using  an 
emulator  routine  which  makes  the  terminal  look  like  a  Tektronix  4014 
to  the  CYBER. 

2.  An  HP9845  is  connected  to  the  HP1000.  Data  can  be  transferred  from 
the  HP1000  to  the  HP9845  and  then  to  the  CYBER.  The  reverse  transfer 
can  also  be  made. 

3.  This  method  is  slow  and  is  not  recommended  for  large  data  files. 

4.  Directions  for  using  the  HP9843  communications  are  given  in  Appendix 

A. 

5.  The  disc  file  on  the  CYBER  MFA  contains  data  in  the  following  format: 

Record  1  contains 

DT  Time  between  samples 

INW  Number  of  data  words 

Records  2,3. . .contain 

INW  data  words 

All  data  are  in  format  (4E20.14).  If  INW  is  not  a  multiple  of  four,  the 
remainder  of  the  last  record  is  filled  with  zeros. 

B.  Magnetic  Tapes 

1.  At  the  present  time,  this  method  is  still  the  most  efficient  and 
least  time-consuming  way  to  transfer  large  amounts  of  data  to  the 
CYBER.  Interactive  usage  on  MFA  greatly  increases  the  clock  time  for 
a  job  which  uses  tapes.  It  is  faster  and  cheaper  to  read  tapes  on 
MFZ.  The  user  has  the  option  of  reading  his  tapes  on  either  MFA  or 
MFZ. 

2.  Directions  for  using  the  procedures  to  read  specific  digital  tape 
formats  and  an  example  are  given  in  Appendix  B. 
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3.  The  disc  file  created  on  the  CYBER  by  the  program  which  reads  the  HP 
tape  contains  data  in  the  following  format: 

Record  (buffer)  1  contains  header  information 

I TYPE  Type  of  File  (0:  A/D  file) 

NAMEF  File  name  on  the  HP1000  tape 

MSSG  Tape  label 

NCHAN  Number  of  channels 

LRECL  Logical  record  length 

NREC  Total  number  of  records 

IREC  Number  of  data  records  (NREC  -  1) 

SRATE  Sampling  rate 

IPBR  Tape  speed  factor 

DT  Time  between  samples 

Records  (buffer)  2  through  NREC  contain  the  data* 

This  format  is  repeated  for  all  the  files  on  the  HP1000  tape. 

4.  To  read  the  disc  file  containing  the  data,  you  must  use  an 
unformatted  READ  for  the  data  records  for  each  file  transferred. 

The  FORTRAN  statements  required  are  as  follows: 

1 0  READ(u )  ITYPE , NAMEF , MSSG , NCHAN , LRECL , NREC , IREC , SRATE , IPBR , DT 

IF( EOF ( u ) . NE . 0 )  GOTO  100 

DO  20  J=l,  IREC 

READ(u)  INW,  ( IDATA( I) , 1=1 , INW) 


20  CONTINUE 


GOTO  10 


100  STOP 

NOTE:  1)  Since  the  current  maximum  logical  record  length  in  the 

programs  which  create  these  tapes  is  8192  (16-bit  words),  the 
number  of  CDC  60-bit  words  required  is  2183  per  record 
(buffer).  The  actual  limit  in  the  CDC  program  to  read  the  tapes 
is  3000. 

2)  At  the  moment,  only  A/D  tapes  (ITYPE=0)  are  converted. 

Other  types  of  tapes  will  be  developed  later.  If  it  becomes 
necessary  to  read  these  tapes  on  the  CYBER,  appropriate  changes 
will  have  to  be  made  in  the  conversion  program.  The  program 
will  print  a  warning  if  ITYPE  t  0,  but  it  will  not  prevent 
execution . 
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3)  u  is  the  unit  number  you  assign  to  read  your  disc  file  in 
your  PROGRAM  statement, 

4)  INW  is  the  number  of  samples  in  each  record. 

5)  Unformatted  input/output  transfers  data  between  main  memory 
and  an  external  storage  device.  Data  are  not  converted  to  any 
format  when  read  or  written  this  way.  However,  the  user  should 
be  aware  that  the  data,  when  written  in  a  formatted  WRITE,  are 
actually  INTEGER  data. 

Once  the  data  have  been  transferred  to  a  CYBER  disc  file,  they  must  be 
sorted  since,  when  the  digitization  is  done,  all  channels  of  interest  are 
sampled  at  one  analog  tape  position  before  proceeding  to  the  next  position. 
This  is  Step  2. 


Step  2.  Sorting  the  data  into  contiguous  sequences. 


The  data  are  obtained  by  sampling  all  channels  from  the  analog  tape  to  be 
digitized  at  a  given  position  before  moving  on  to  the  next  position,  as 
illustrated  in  Figure  2. 


TAPE  1 
POSITION 
N 

Figure  2. 


jwMBHI 

fl- 

f-B 


e 


I  tape  I 

POSITION 

N+1 

Data  Sampling  Format 


data  oharma!  1 


data  oharmal  2 


date  ohormat  J 


Since  the  analog  tape  is  actually  moving  at  a  constant  speed,  there  is  a 
slight  time  skew  from  channel  to  channel,  but  this  skewness  is  normally 
negligible  compared  to  the  sampling  rate. 

A  computer  program  has  been  written  for  MFZ  to  sort  the  data  into 
contiguous  sequences,  plot  the  sequences,  and  write  a  permanent  file  for  the 
next  step.  (See  Appendix  C.)  Since  all  of  the  data  for  one  channel  are  not  in 
memory  at  any  one  time,  several  sequences  are  formed  and  stored  on  temporary 
files,  one  file  for  each  channel,  until  the  end  of  the  digitizing  pass.  Each 
temporary  file  is  then  written,  in  order,  to  the  permanent  file  in  such  a  way 
that  all  of  the  sequences  for  the  first  channel  are  written  first,  all  of  the 
ones  for  the  second  channel  are  next,  and  so  forth.  This  procedure  is 
illustrated  in  Figure  3. 
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INPUT  FILE 


DATA  :  X  .  .  . 

J.  1.  k 

j  -  1.2 . J  (RECORD  NUMBER) 

l  «  1.2 . 1  (ITEM  NUMBER  WITHIN  RECORD) 

k  -  1,2, ...,K  (DATA  CHANNEL  NUMBER) 


TEMPORARY  FILES 


OUTPUT  FILE 


NOTE:  Temporary  files  are  necessary  because  the  data  are 
read  one  record  at  a  time  and  the  quantity  of  data  may  be 
very  large. 


Figure  3.  Schematic  of  the  Procedure  To  Sort  the  Data 
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In  the  example  shown  in  Figure  4,  the  calibration  steps  and  the  event  are 
digitized  in  the  same  pass.  This  method  creates  an  excess  of  data  points.  An 
editing  procedure  to  save  only  the  pertinent  data  is  discussed  next  in  Step  3. 

Step  3.  Editing  the  Data  File 

A  typical  channel  of  data  can  be  divided  into  seven  sections,  as  shown  in 
Figure  3.  Not  all  channels  contain  every  section.  Sections  1,  3,  5,  and  7 
contain  data  which  can  be  deleted;  however,  the  duration  of  Section  5  must  be 
accounted  for  if  time  is  counted  from  the  fiducial  mark  in  Section  4  and  is 
stored  in  the  variable  TSTART  on  the  output  file.  Using  this  rationale  on  the 
example  in  Figure  4,  sequences  1,  2,  13,  17,  18,..., 23  can  be  deleted.  Now, 
sequences  3  through  12  make  up  section  2;  sequence  14  makes  up  section  4;  and, 
sequences  13  and  16  make  up  section  6.  There  is  no  section  5  in  this 
example.  Next,  we  sample  section  2  at  equally  spaced  intervals.  Thus,  using 
this  method,  we  have  vastly  reduced  the  number  of  data  points. 

A  computer  program  has  been  written  for  MFZ  to  edit  the  data,  plot  the 
results,  and  write  a  permanent  file  to  be  used  in  Step  4,  where  the  conversion 
to  engineering  units  is  done.  The  input  parameters  required  for  the  program 
are:  1)  the  sequence  numbers  where  sections  2  through  7  start,  2)  the  sampling 
interval  for  the  calibration  steps  (for  instance,  if  one  sample  is  taken  for 
every  fifty  input  data  points  -  the  interval  equals  fifty),  3)  the  sampling 
interval  for  the  event,  and  4)  the  option  to  print  a  message.  For  the  example 
in  Figure  5,  the  input,  format  915,  would  be: 

3  13  14  15  15  17  50  1  0 

Following  this  process,  the  data  points  for  a  given  channel  can  finally  be 
merged  into  one  sequence.  A  computer  listing  of  Step  3  is  in  Appendix  D.  A 
plot  of  the  results  thus  far  is  shown  in  Figure  6. 

Step  4.  Converting  to  Engineering  Units 

A  computer  program^  was  written  several  years  ago  to  convert  digitized  data 
into  engineering  units.  The  program  was  modified  to  accept  the  data  produced 
in  Step  3.  The  mechanics  of  handling  the  data  are  different,  but  the  actual 
conversion  is  the  same.  The  equipment  producing  the  data  and  recording  them 
on  analog  tape  has  not  changed;  just  the  digitizing  equipment  has  changed. 

This  FORTRAN  IV  program  directs  the  digital  computer  to  read  the  information 
from  the  computer  file,  perform  the  many  necessary  calculations,  and  provide 
the  reduced  data  in  tabular  and/or  plotted  form  as  output  of  the  computer. 

This  program,  which  is  in  modular  form,  has  many  options  which  are  chosen 
according  to  the  particular  input  data  and  the  desired  output.  The  options 
available  to  date  in  the  program  are:  (1)  calibration,  (2)  location  of  a  step 
function,  (3)  conversion  of  the  dependent  variable  to  engineering  units,  (4) 
integration  of  the  dependent  variable,  (5)  plotting  of  the  output  data,  and 
(6)  tabulation  of  the  output  data. 


^C.L.  Henry,  /?.L.  Martz,  E.M .  Wineholt ,  "An  Improved  Procedure  For  The 
Reduction  of  Interior  Ballistic  Data  Recorded  on  Analog  Tape BRL  MR 
2374  ,  April  1974  (AD  919924L). 
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Figure  5.  Typical  Sections  of  a  Data  Channel 


Figure  6.  An  Example  of  the  Data  Remaining  After  Editing 

Basically,  the  conversion  is  accomplished  by  fitting  a  second-degree 
polynomial  by  the  method  of  least  squares  to  the  amplitudes  in  counts  of  the 
calibration  steps  which  have  known  values  in  engineering  units.  Thus,  the 
relationship  between  engineering  units  and  counts  is  determined  and  the 
conversion  can  be  accomplished.  The  time  for  each  data  point  relative  to  the 
fiducial  mark  is  also  calculated.  The  option  to  integrate  once  or  twice  is 
provided  as  well  as  the  option  to  plot  or  not.  A  listing  of  the  program  and  a 
sample  input  and  output  are  in  Appendix  E. 
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Various  parameters  are  necessary  input  data  for  the  program,  as 
illustrated  in  Figure  7.  They  are: 


1. 

NX 

The  number  of  samples  to  be  averaged  for  each 
calibration  step. 

2. 

NY(i),i=l,...,NS 

The  sample  indices  marking  the  position  on  each 

(NS  <  6) 

calibration  step  to  start  processing  data. 

3. 

ITZ 

The  number  of  samples  to  skip  after  sampling  the  last 
calibration  step  before  starting  to  search  for  the 
fiducial  mark, 

4. 

IX 

The  number  of  samples  to  skip  after  the  fiducial  mark 
to  reach  the  data  which  are  to  be  converted  to 
engineering  units, 

5. 

IBSE 

The  number  of  samples  to  skip  after  the  fiducial  mark 
before  starting  to  sample  the  baseline, 

6. 

B 

The  calibration  constant  for  the  gage  in  engineering 
units  per  calibration  step,  and 

7. 

DLTM 

A  time  adjustment  to  be  subtracted  from  the  value  at 
the  fiducial  mark,  TSTART,  which  comes  from  the  data 
file. 

Some  channels  of  data  may  not  have  a  fiducial  mark  so  the  program  has  the 
option  to  skip  that  part.  Time  is  then  counted  from  the  first  sample. 

These  parameters  and  other  control  variables  either  have  default  values 
which  can  be  changed  using  NAMELIST  or  are  read  in  as  card  images.  The 
details  and  an  example  are  provided  in  Appendix  E. 


— »1NXK — 


Figure  7.  Graphic  Explanation  of  Input  Parameters 


III.  CONCLUSIONS 


To  convert  experimental  data  digitized  by  IBD's  HP1000  minicomputer  system 
into  the  same  computer  file  format  required  by  existing  computer  programs,  a 
four-step  procedure  has  been  devised.  Separate  steps  are  necessary  because 
engineering  judgments  must  be  made  before  further  processing  is  possible. 

These  steps  are: 


1. 

Reading  the  data  into  a  computer  file, 

2. 

Sorting  the  data  so  that  all  the  data  samples  for 
contiguous  in  a  file. 

one 

channel  are 

3. 

Editing  the  data  so  that  unnecessary  samples  are 
further,  and 

not 

processed 

4. 

Converting  the  data  into  engineering  units  and  creating  an  output 
file  with  the  same  format  as  previously  generated  in  Reference  1. 

This  output  file  now  contains  the  data  in  the  format  required  for  input  into 
numerous  experimental  data  analysis  programs . 


IT 


APPENDIX  A 

INPUT  TO  CYBER  MFA  VIA  MODEM  LINK  ON  HP9845 
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APPENDIX  A 


Preparation . 

a.  Turn  on  HP9845  and  7906  disc  drive. 

b.  When  disc  drive  is  ready, 

LOAD  "ANPACK:D12" , 10 


IEXECUTEI 


c.  When  menu  appears, 

K3 1  (Data  Communications) 


d.  Choose  data  link. 


for  HP1000 


2  ICONT  I  for  CDC 


e.  Specify  HP9845  file. 

When  "terminal  ready"  message  appears,  indicating  the  terminal 
emulator  program  is  loaded,  then 

|K13|  (1)  To  get  into  edit  mode 

(2)  It  is  not  necessary  to  change  any  info  on  the  first  line. 

If  you  are  going  to  the  CDC  and  you  do  not  want  your 
password  displayed  on  screen,  you  may  want  to  change  ECHO 
to  ON  by  using  |  >  |  to  space  over  and  use  STEP  key  to 

change  to  'ON'. 

STOReI  -  To  store  data  communication  information  line 

(3)  This  line  needs  to  be  changed  to  the  name  of  file  and 
its  size  on  the  HP9845. 

Change  TEST:TI5  to  filename : C 1 2 

Change  SIZE=010  to  #  of  records  you  want 

STORE  -  To  store  file  info  line 

(4)  Will  get  a  message  that  edit  mode  has  been  exited. 

(5)  If  this  is  a  new  file,  then  it  must  be  created.  If  it 
is  not,  skip  this  step. 

ISHIFTlnon  To  create  file  on  HP9845 

Answer  "Y"  to  creation  question.  Message  will  appear  on 
screen  when  creation  is  completed. 
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2. 


3. 


f.  CONT  [Sends  prompt  to  CDC  or  HP1000.  Wait  for  LOGIN  messages. 


From  CDC  to  HP9845. 

a.  LOGIN  to  CDC 

b.  Repeat  I.E.(l)  through  I.E.(4),  changing  first  edit  line  back  to 
ECHO  OFF  if  you  turned  it  ON. 


/ GET, A=Pfn . 
/ASSIGN, TT,B. 
/COPYSBF,A,B. 

EH 


CONT 


SHIFT 


(CK 

No  fCR)  ! 
to  turn  record  on 


CONT 


When  finished  recording, 

|K15l  to  turn  record  off 


SHIFT 


Log  off  CDC 
/BYE 

l~K4~l 


SHIFT 


to  disconnect  from  CDC; 


now  back  under  control  of  HP9845. 


From  HP9845  to  CDC 

a.  LOGIN  to  CDC 

b.  Repeat  I.E.(l)  through  I.E.(4),  changing  first  edit  line  back 
to  ECHO  OFF  if  you  turned  it  ON. 

c.  Terminal  definition  defaults  are  set  to  emulate  a  Tektronix 
4014  with  a  page  width  of  80  columns.  If  you  need  to  set  the 
terminal  definition  (TC  parameter)  to  another  terminal  or  the 
page  width  to  a  different  size,  then  enter 

/TRMDEF.TC*  ,PW« 

See  the  NOS  IAF  Manual  for  parameter  values.  For  most  data  transfers 
from  the  HP9845 ,  the  default  parameters  will  be  good  and  this  step  can 
be  skipped. 

d.  Data  transfer  procedure: 

NEW, LFN1 
TEXT 

'Enter  Text  Mode'  message  returned  by  CDC 
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4. 


SHIFT 


ms 


"N11  to  handshake  message  from  the  HP9845 

Data  being  transferred  will  be  displayed  on  screen. 

When  all  data  has  been  transferred, 

m  then 


tONTROLl 


CONT 


This  terminates  the  input  in  TEXT  mode  and  gets  you  back  to  READY 
mode  in  IAF . 


e.  Check  the  data  transfer: 
/REWIND, LFN1 

/LIST, F=LFN1 

f.  To  save  the  file: 

/ SAVE, LFN1=PFN/PW= 


,M= 


,  CT= 


g.  Disconnect  procedure: 
/BYE 
|SHIFT 


fK4l 


From  HP9845  to  HP1000. 

a.  Log  in  to  HP1000. 

b.  Repeat  steps  I.E.(l)  to  I.E.(4),  changing  file  info  line  to 
correct  file  name  and  file  size. 


:ST,1,NAMR  -  No  (Ci 
|SHIFT|  (¥171  Upload 

"N"  to  handshake  question 


CONT 


CONT 


d.  When  finished  transferring  to  HP1000, 

DO  NOT  |SHIFT|  I K14]  as  this  will  cause  great  problems. 

e.  Be  careful  not  to  hit 


SHIFT 


f  K1 5 1  as  this  will  write  all  kinds  of  error  messages 


over  your  file  on  the  HP9845. 
Log  off  HP1000 
:EX,SP 
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g.  (shift 


dD 


to  disconnect. 


From  HP1000  to  HP9845. 

a.  Log  in  to  HP1000. 

b.  Repeat  steps  I.E.(l)  to  I.E.(4),  changing  file  info  line  to 
correct  file  name  and  file  size. 


C.  ST,  NAMR, 1  No  @ 

I K14 1  to  turn  record  on 
to  start  recording 


SHIFT 


CONT 


d.  When  finished  recording, 

| K14 |  to  turn  record  off 


SHIFT 


e.  Log  off  HP 1000 
:EX, SP 

tm 


f.  SHIFT 


to  disconnect 


CONT 
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APPENDIX  B 

INPUT  TO  CYBER  VIA  MAGNETIC  TAPE 
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APPENDIX  B 


1.  MFA 

Usern . 

USER, username  password • 

CHARGE, account , pro ject  number, 

GET, HPTAPEB. 

F ILE , TAPE7 , CM=NO . 

BEGIN , GET , TAPE , LF-TAPE  7 , VSN-vsn , DEN= 1600, TK=9 , LABELED-NO , S . 

where  vsn=tape  library  number 
HPTAPEB . 

REWIND, TAPE 10, 

SAVE , TAPE lO^pfn • 

where  pfn  =  name  of  disc  file  you  are  creating. 

Note  that  input  tape  from  the  HP1000  must  be  on  unit  7  and  your  output  disc 
file  is  on  unit  10. 

2.  MFZ 


Usern, STMFZ,P4, NT 1 . 

ACCOUNT, account . 

BEGIN , READZ , HPTAPEP , VSN=vsn , PF=pf n , ID=user id . 

where  vsn  =  tape  library  number, 

pfn  ~  name  of  disc  file  you  are  creating,  and 
userid  =  your  ID  on  MFZ. 

It  is  strongly  recommended  that  you  use  this  procedure  and  copy  your  file  from 
MFZ  to  MFA  if  you  need  the  data  on  MFA. 

3.  Listing  of  MFA  file  ADTAPE1/UN=B00TS: 

BOOTS, STMFZ,P6,NT1 . 

ACCOUNT, PDxxx. 

BEGIN , READZ , HPTAPEP , VSN=69  7  8 , PF=L0LA2 , ID=B00TS . 
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APPENDIX  C 

IMPLEMENTATION  OF  STEP  2 
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APPENDIX  C.l.a 

Listing  of  MFA  File  ADTAPE2/UN=B00TS 
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APPENDIX  C. l.a 


Listing  of  MFA  File  ADTAPE2/UN=B00TS 

a.  Job  Control  Language 

BOOTS , STMFZ , P6 ,T100 . 

ACCOUNT ,PDxxx. 

BEGIN , ATTACH , PLOTLIB . 

ATTACH , TAPE 1 , L0LA2 , ID=BOOTS . 

REQUEST , TAPE3 ,  *  PF . 

REQUEST, TAPE2,*PF. 

REQUEST, TAPE13,*PF. 

FTN,R=iO. 

LGO. 

EXIT, U . 

BEGIN, PLOT, CALCOMP, TAPE 13 . 

CATALOG , TAPE3 , LOLAD , ID=BOOTS . 

CATALOG , TAPE 1 3 , PLTT , ID=BOOTS . 
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APPENDIX  C.l.b 
Program 
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1 


6 

It 

15 


2# 

25 

31 

35 


PRO GRAN  RDNTD  (INPUT*  OUTPUT*  TAPES*  I  NPUT*TAPES*  OUTPUT*  TAPE  I  *TAPE3* 

*  TAPE13*TAPEH,TAPE1S,TAPE16*TAPE1?*TAPE1«*TAPEI#* TAPES* 
t  TAPE2#* TAPE21* TAPE22*  TAPE23.TAPE84* TAPESS*  TAPES# ) 

CONNON  I  STOP  (2500  )1C0UHT(85M)*NCH*NIUC*NU*E*lR01CHl(25##)eIE*ISU* 
t  CH8(SS##)*CH3(&5##)*CH4(25##)*CHS(3S##)*CHi(85##)iCH?(35##)<NCV* 

I  CH8(85##)*CH9(S5##  )*CH1#(2S##2#CH11(8S##)*CH12(8S##)*CAL(8S##) 
t*NCAL*NCHA*N$SG(8  )*  IPL*IDIN#NANEF*DT*  JJ*NUS 
IDIN-SS## 

IPL-1 

IRD-# 

5#  UR  I  TEC  6*  IMS ) 

CALL  RDDTA(l) 

DO  6#  IU-14,26 
REWIND  IU 
6#  CONTINUE 
IRD-IRD+1 
NBR-NCH*(IDI N/NU) 

NCV-NBLK/N8R 

IF ( NILX-NCYtNBR . NE . • )  NCV*NCY*1 

ISU-i 

IE-1 

NWS-# 

NCAL-# 

NCHA-# 

DO  5##  I-1*NCV 

IF  ( I  * EQ .HCY .  AND .  I  A»S ( N1UC-NCY8N8R ) . HE •  #  )NiR- 1 AiS ( N8LK-  ( NC V- 1  ) 8NBR  ) 
DO  1##  N-1*NBR 
CALL  RDDTA(2) 

CALL  DENULT(N) 

1##  CONTINUE 
JC-K-l 

CALL  WRTDTA 
5##  CONTINUE 

URITEC6, 11## )  NCAL*NCHA 
CALL  URTTPE 
GO  TO  5# 

1###  FORNATC 1H1 ) 

11##  FORNATC  *  NO.  OF  CAL  BLKSM1#*SX*  'NO.  OF  DATA  1USM1#) 

END 


37 


1 


6 

It 

15 

30 

85 

30 

35 

40 

45 

55 

55 

60 


SUBROUTINE  DENULTCICVL) 

COHHON  IST0RCB6OO)»C0UNTC8S00>»NCH,NBLK,NU,K,  IRD.CH1  C8S00 ).IK,  ISU 
*,CH8C3SO0>,CH3C8500>»CH4C3S00)»CHS<86OO)»CHSC8S00)»CH7C8S00)»NCV, 
8  CH8C3SO0>,CH9<aS00>.CH10<3S0OKCNllC8SO0>tCH13<a60OI.CAI.CB*60) 
8,NCAl,NCHA,H$SG(8 ), IPL,IDIA,NAN6F,DT, JJ»NUS 
IFCICVL.EQ.DJJ-O 
IFCRJ.EQ..0>JJ-0 
IFCJJ.E0.0)HS-1 
IFCICYL.EQ.l  )K«1 
HUS-NUS+NU 
RJ»HODCNU$»NCH) 

IF(JJ.CT.B)OOTOilS 
90  DO  100  H-NS,NU,NCH 

IF(RJ.GT.O. .AND.N+HCH.GT.NU >  00  TO  110 

IFCISO.EQ.l)  IK-K 

XSU-8 

IF(NCH.CE.l)  CHICKWSTORCH) 

IFCNCH.GE.3)  CH3CK>-IST0RCH*1 ) 

IF(NCH.GE.3)  CH3(K)*IST0R<n+3 ) 

IF ( NCH .  GE .  4 )  CH4CK)-IST0RCN*3> 

IF CNCH.GE.5)  CH5CK)"IST0RCN*4) 

IFCNCH.GE.6)  CH6(K)-I$T0R(H45) 

IF(NCH.GE.T)  CH7(IC)«IST0R(H46) 

IFCNCH.GE.8)  CH8CK)-IST0RCN47) 

IF(NCH.GE.S)  CH9(K)-XST0R(N48) 

XF(NCH.GE.IO)  CH10CK)«IST0RCH+9> 

IFCNCH.GE.il )  CH11CK>«XST0RCN+10) 

XFCNCH.GE.18)  CH13CK)-IST0RCH+11 ) 

GO  TO  60 

50  CAKE)*  1ST 0R(H> 

60  K-K41 

IFCK.EO.il .AND. ISU.EQ.l )  URXTEC6, 1000 )  IRD,  CCALC J ), J-l, 10) 
IFCK.EQ.il .AND. ISU. EQ. 3)  URXTEC6,1000)  IRD,  <CH1CJ),J-1.10) 

100  CONTINUE 
RETURN 

110  JJ*NW-H+1 
NN-NU-JJ+1 
GO  TO  130 
115  NH—JJ+l 
JS-JJ 


JI-JJ41 

JJ'NCH 

NS-NCH-JS+1 

GO  TO  C 130,131,132, 133,184,135,136,137, 138,189,130, 131 ), JI 

130  XFC  JJ.GE.l )  CHICKWSTORCHN) 

131  IFC JJ.GE.3)  CH2CK >*I$T0R<HN+1 ) 

133  IFC  JJ.GE.3)  CH3CK)-IST0R(m«3) 

133  IFC JJ.GE.4)  CH4(K)«IST0RCnn+3> 

134  IFCJJ.GC.6)  CHSCK)-ISTOR(NNM) 

135  IFC JJ.GE.6)  CH6CK )aIST0RCNN+5 ) 

136  IFCJJ.GE.7)  CH7CK)-IST0RCNN46) 

137  IFCJJ.GE.8)  CHS<K)«IST0RCnn+?) 

138  IFCJJ.GE.9)  CH9<K)sIST0RCHN+8) 

139  IFC JJ.GE. 10)  CH10CK)-IST0RCWH9> 

130  IFCJJ.GE.il)  CHUCO-ISTORCWHIO) 

131  IFC JJ.GE. 13)  CM13CK>-IST0RCRHM1> 

IFC JJ.EQ.NCH)KSK+1 


IFC JJ.EQ.NCH)GO  TO  90 
RETURN 

1000  F0RNATC1H  ,I6,10FS.0> 
END 
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SUBROUTINE  PLAT<  IPLOT,DATA,IRDl,  INU ) 

COHNON  IST0R(2$99 ),  COUNT  <2599  ),NCH,NBLK,  NU,K,  I RD,  CHI  (2599  ),IIC,J^y 
*TCM8(2599)#CH3(2599)#CH4(B599>#CH5(2S90>,CH6<BS09>,CH7<8S99>,NCY# 
*CH8 ( 2599  >  ^  CHS  ( 2599  ),CH19( 2599),  CHI  l  (2599),  CHIB  (8599),  CAL  (8599) 
B.NCAL/NCHAjNSSC (B ),  I^JrilDINfNAHEF,  DT,JJ,NW$ 
fclHENSION  LABEL (4  ),  DATA(2S99  >,  ITITLEO ) 

IF ( IPLOT.GT . 1 )  CO  TO  199 


IPLOT-2 
XB-1.75 
VI*  .4 
XAX*5. 

VAX* 4. 

XPAGE*7* 

YPAGE  *4*75 
FACT-1. 

IUNIT-13 

LABEL ( 1 )-L9H  BOOTS  399 
LABEL<2)*19H  6121 
LABEL ( 3  ) * 19HRAU  DATA 
LABCL<4 )-19HL0LA2 
NODE-1 

CALL  PLTBEG(XPAGE,YPAGE,FACT,IUNIT,LAB£L> 

XNI  *9. 

XNA-IDIN 

YN1-9. 

YNA-5999. 

DX-599. 

DV*  1999. 

XS-(XRA-XHI)/XAX 

VS-(VHA-VHI)/VAX 

199  CALL  PLTSCA(XB,YB,XNI,YNI,X$,VS) 

CALL  PLTUND<XNI.XHA,YHI,YNA) 

CALL  PLTAXSCDX,  DY,XNI,XRA, YNI, YHA,4  ) 

CALL  LA9€LA(DX,BV#XHI#XHA/VflI,VHA/l.#l.  ) 
ENCODE <39# 1999# I TITLE)  IRD1,INU 
1999  F0RNAT(3H  ID,I19,8H  PL0T#I8#1H>) 

TX-XNI 

TY-YNA+.95BYS 

CHT-.l 


CALL  PLTSYH(CHT, ITITLE (1  ),9. #TX,TY) 

TIC -DXBDTS 1999. 

TINE-DTSFLOAT ( NM )tl999 . 

TY-YNA+.20IYS 

ENCODE ( 12#2999#  ITITLE  )NANEF 
2999  FORNATdH  ,A19,  1H>  ) 

CALL  PLTSYH<CHT, ITITLE, 9. ,TX,TY) 

TX-XNI-XStl.SS 

TY-YNA-YSI.3 

ITITLE(1  )«19HTIHE,  RS> 

CALL  PLTSYmCHT,  ITITLE, 9. ,TX,TY) 
TY*YHA-YSS.4S 

ENCODE ( 16, 3999/ITITLE )TIHE 
3999  FORNAT ( 8H  /FRAK-.F7.3, 1H>  ) 

CALL  PLTSYH(CHT, ITITLE, 9. ,TX,TY) 
TY-YHA-YSt.6 

ENCODE < 16,4999, ITITLE  )TIC 
4999  F0RNAT(8H  /TIC*  ,F7.3,1H>) 

CALL  PLTSYH(CHT,ITITLE,9.,TX*TV) 

CALL  PLTD TS ( NODE , 9 , COUNT ( 1 ) , DATA ( 1 ) , NU , 9 ) 

CALL  PLTPGE 

RETURN 

END 
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COHNON*ISTOR?2SOO),COUNT(2SOO),NCH,N8LK,NW,K#IRD»CH1 (2500  )> IK,I$W 

*  CH2(SSO0>,CH3(25OO),CH4(2SOO),CHS(2SO0 ), CH6( 2500 ),CH7( 2500 ),NCV, 

*  CH8(25O0>,CH9(25OO),CHtO(86OO),CHll(8SO0 >,  CHI 8 (8500),  COL  (8500 ) 
*.NCOl,NCHA,NSSG(8>,IPt.,IDin,NOnEF,DT,JJ.NyS 

fFdSW.EG.S)  00  TO  40 
WRITE (14)  K,(CAl(J),J-l.K) 

NCAL-NCAl+1 

RETURN 

40  IFdK.EQ.l.AND.  ISW.E0.2)  00  TO  SO 

IF(KK~NE.O>WRITE(  14 )  KK,(CAKJ),J-1.KK> 

IF(KK.NE.O)NCOl-NCRld 


IJ-K-KK 

IF(NCH.CE.l)  WRITEdS)  I J, (CHI ( J ), J-IK 
IF(NCH.CE.2)  URITE( 16 )  I J, (CH2( J ), J-IK 
IF(NCH.GE.3)  URITE( 17)  IJ. (CH3( J),J-IK 
IF(NCH.CE.4)  URITE( 18 )  IJ. (CH4( J ), J-IK 
IF(NCH.CE.S)  WRITEdS)  IJ.(CHS( J),J-IK 
IF(NCH.CE.6)  URITEOO)  IJ,  <CH6<  J >,  J-IK 
IF(NCH.GE.7)  URITEC21 )  IJ,<CH7( J),J-IK 
IF(NCH.0E.8)  URITEC22)  I J, (CH8( J ), J-IK 
IF(NCH.GE.S)  WRITE (23)  I J, (CHS( J ), J-IK 
IF( NCH.GE. 18 )  WRITE (24)  IJ, (CH10( J ), J- 
IF(NCH.OE.n)  WRITE(2S)  IJ, (CH11( J), J- 
IF(NCH.GE.12>  WRITE (26)  I J, (CH12( J ), J- 
IK-1 

NCHO-NCHftd 

RETURN 


K) 

K) 

K) 

K) 

K  > 

K) 

K) 

K) 

K) 

K,K ) 
K,K) 
K,K) 


50  IF(NCH.GE.l)  WRITEdS)  K,  (CHI  ( J ),  J-1,K ) 
IF(NCH.0E.2>  WRITE( 16)  K, (CH2( J ), J-1,K ) 
IF(NCH.GE.3>  URITE( 17  )  K, (CH3( J ), J-1,K  ) 
IF(NCH.GE.4)  WRITE (18)  K, (CH4( J ), J-1,K  ) 
IF(NCH.GE.S)  URITE(19)  K, (CH5(J ), J-1,K  ) 
IF(NCH.GE.6>  URITE(20)  K, (CH6( J ), J-l.K  ) 
IF(NCH.GE.7)  WRITE(21 >  K, (CH7( J), J-1,K  ) 
IF(NCH.GE.8)  WRITE (22 )  K, (CH8( J), J-1,K  ) 
IF(NCH.GE.S)  WRITE (23)  K, (CH9( J ), J-1,K ) 
IF(NCH.GE.IO)  WRITE (24 )  K, (CH10( J ), J-l ,K ) 
IF (NCH.GE. 11 )  WRITE(2S)  K,(CH11(J),J-1,K) 
IF(NCH.CE. 12 )  WRITE(26 )  K, (CH12( J), J-1,K ) 
WRITE (6,1 >NCV,NCHA,K, ( CHI (J), J-l, 10) 

1  F0RNAT(3I10,10F7.0) 

HCHO-NCHR+1 

RETURN 

END 
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1 


5 

It 

IS 

85 

3# 

35 


45 


subroutine  urttpe  _ 

COHNON  I STOR< 2500 ). COUNT (8500 ),NCH,NOLK,NU,K, IRD.CM1  (8600).IK.  18U, 

*  CH8<8SOO>,CH3(2SOO).CH4<8500>,CHS<8SOO>.CH«<8SOO).CH7<8500>.HCY, 

*  CH8(8SOO>,CH»C8SOO),CH10<8SOO),CNll<aSOO),CH18<aSOO>,CAl(8SOO) 

*  ,NCAL,NCHA,HSSC(0),IPL, IDIH,NAHEF,DT, JJ.NU6 
DIHENSION  DATA<8500) 

DO  5  1*1, ID1H 
COUNT ( I )*I 
5  CONTINUE 
ITPE-14 
IPL*1 

IF ( NCAL .  EO .  0 )00  TO  300 

URITE(2)NAHEF 

IRD1-IRDMOOOO 

UNITE (8)  IRDl.DT.NCH.NCAL 

REWIND  ITPE 

DO  800  H*1,HCAL 

READ(ITPE)  NU,(DATA<J),J*1,NU) 

IF(N.EQ.l)  URITE<6,1000)  IRD1,HU,<DATA<J),J-1,10> 

IFCN.E0.NCAU8)  URITE<6,1000>  IRDl,NU,<DAtA<J). J-1,10) 
IF(N.EO.NCAL)  URITE(6, 1000)  IRDl.NU.  CDATAU),  J-1,10) 

IF(N.EQ.1)CALL  PL AT ( IPL , DATA. IRD1 , N ) 

URITEC2)  HU, ( DATA (J),J*1,NU) 

800  CONTINUE 
300  ITPE*ITPEfl 
URITEt  3 1NANEF 
WRITEC3)  IRD.DT.NCH.NCHA 
DO  500  N-l.NCH 
IRD1*IRD+100*H 
REUIND  ITPE 
DO  400  I-l.NCHA 

READ(ITPE)  NU,<DATA<J),J*1,HU) 

IF ( I .EO . 1 >URITE< 5, 1000 >ITPE»NCHA 
IFd.EO.l  )URITE(6,1000)NCH 

IF ( I .EO. 1 )  URITE(6,1000)  IRDl.NU, (DATA< J),J*1,10) 

IF( I.EQ.NCHA/2)  URITEC6.1000 )  IRDl.NU, <DATA<J), J-1,10) 
IF(I.EQ.NCHA)  URITE(6,1000)  IRDl.NU. (DATA(J).J-l, 10) 

IFCN.EO.l >CALL  PLAT  (IPL.DATA.IRDl.l) 

WRITE(3>  HU, (DATA(J), J*1,NU) 

400  CONTINUE 
ITPE-ITPEM 
600  CONTINUE 
RETURN 

FORNATdN  ,2110, 10F8.0) 

END 
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APPENDIX  C.2 
Sample  Output 
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FILE  NANE  L2R47 
NO.  OF  RECORDS  • 
TINE  INCRENENT,SEC 


NO.  OF  CHANNELS  ■ 
288  LENGTH  OF  A  RECORD 

12S88E-94 


4 


1024 


2188. 

2165. 

2162. 

2163 

2162.  2168. 

2158. 

2158. 

2163. 

2167. 

25 

e 

2848 

2168. 

2165.  2162.  2163 

.  2162. 

2168. 

2158. 

2158. 

2163. 

2167 

2142. 

2143. 

2144. 

2148 

2148.  2143. 

2151. 

2155. 

2158. 

2143. 

25 

1 

2848 

2142. 

2143.  2144.  2148 

.  2148. 

2143. 

2151. 

2155. 

2158. 

2143 

2121. 

2123. 

2124. 

2124 

2119.  2118. 

2118. 

2123. 

2129. 

2127. 

25 

2 

2848 

2121. 

2123.  2124.  2124 

.  2119. 

2118. 

2118. 

2123. 

2129. 

2127 

2880. 

2894. 

2896. 

2895 

2898  .  2891. 

2899. 

2188. 

2897. 

2892. 

25 

3 

2848 

2889. 

2894.  2896.  2895 

.  2898. 

2891. 

2899. 

2189. 

2897. 

2892 

2358. 

2359. 

2368. 

2363 

2363.  2364. 

2368. 

2371. 

2367. 

2365. 

25 

4 

2848 

2358. 

2359.  2368.  2363 

.  2363. 

2364. 

2368. 

2371. 

2367. 

2365 

2368. 

2363. 

2361. 

2368 

2366.  2374. 

2375. 

2369. 

2362. 

2369. 

25 

5 

2848 

2368. 

2363  .  2361.  2368 

.  2366. 

2374. 

2375. 

2369. 

2362. 

2368 

2645. 

2639. 

2636. 

2635 

2634.  2634. 

2637. 

2643. 

2646. 

2644. 

25 

6 

2848 

2645. 

2639  .  2636  .  2635 

.  2634. 

2634. 

2637. 

2643. 

2646. 

2644 

2636. 

2646. 

2658. 

2646 

2639.  2632. 

2625. 

2624. 

2628. 

2634. 

25 

7 

2848 

2636. 

2646  .  2658  .  2646 

.  2639. 

2632. 

2625. 

2624. 

2628. 

2634 

2988. 

2982. 

2983. 

2986 

2912.  2919. 

2922. 

2917. 

2913. 

2911. 

25 

8 

2848 

2988. 

2982.  2983.  2986 

.  2912. 

2919. 

2922. 

2917. 

2913. 

2911 

2916. 

2916. 

2911. 

CWW9 

2988.  2898. 

2984. 

2912. 

2915. 

2917. 

25 

9 

2848 

2916. 

2916.  2911.  2985 

.  2988. 

2898. 

2984. 

2912. 

2915. 

2917 

3193. 

3193. 

3198. 

3187 

3184.  3185. 

3186. 

3186. 

3186. 

3188. 

25 

18 

2848 

3193. 

3193.  3198.  3187 

.  3184. 

3185. 

3186. 

3186. 

3186. 

3188 

3184. 

3188. 

3198. 

3191 

3194.  3195. 

3192. 

3186. 

3181. 

3181. 

25 

11 

2848 

3184. 

3188.  3198.  3191 

.  3194. 

3195. 

3192. 

3186. 

3181. 

3181 

3467. 

3466. 

3468. 

3469 

3468.  3465. 

3463. 

3462. 

3468. 

3461. 

25 

12 

2848 

3467. 

3466.  3468.  3469 

.  3468. 

3465. 

3463. 

3462. 

3469. 

3461 

3464. 

3461. 

3463. 

3467 

3478  .  3469. 

3469. 

3469. 

3468. 

3463. 

25 

13 

2848 

3464. 

3461.  3463.  3467 

.  3478. 

3469. 

3469. 

3469. 

3468. 

3463 

2894. 

2692. 

2891. 

2897 

2897  .  2897. 

2894. 

2894. 

2899. 

2896. 

2896 

25 

14 

2848 

2894. 

2892.  2891.  2897 

.  2897. 

2897. 

2894. 

2894. 

2899. 

2889. 

2682. 

2879. 

2885 

2886.  2885. 

2882. 

2882. 

2889. 

2893. 

2893 

25 

15 

2848 

2889. 

2882  .  2879  .  2885 

.  2886. 

2885. 

2882. 

2882. 

2889. 

2882. 

2683. 

2888. 

2887 

2888.  2877. 

2878. 

2885. 

2884. 

2877. 

2877 

25 

16 

2848 

2882. 

2883  .  2888  .  2887 

.  2888. 

2877. 

2878. 

2885. 

2884. 

2366. 

2366. 

2378. 

2371 

2373.  2375. 

2374. 

2378. 

2373. 

2378. 

2378 

25 

17 

2848 

2366. 

2366.  2378  .  2371 

.  2373. 

2375. 

2374. 

2378. 

2373. 

1 

2367. 

2368. 

2365. 

2363 

2367.  2369. 

2369. 

2367. 

2369. 

2372. 

2369. 

2372 

25 

18 

2848 

2367. 

2368  .  2365  .  2363 

.  2367. 

2369. 

2369. 

2367. 

1 

2299. 

2368. 

2319. 

2338 

2338.  2344. 

2352. 

2358. 

2362. 

2365. 

2365 

25 

19 

2848 

2299. 

2388.  2319.  2338 

.  2338. 

2344. 

2352. 

2358. 

2362. 

1 

2376. 

2377. 

2377. 

2371 

2367  .  2369. 

2373. 

2376. 

2378. 

2377. 

2377 

25 

26 

2848 

2376. 

2377  .  2377  .  2371 

.  2367. 

2369. 

2373. 

2376. 

2378. 

1 

2385. 

238S. 

2385. 

2382 

2374  .  2368. 

2369. 

2378. 

2383. 

2375. 

2375 

25 

21 

2848 

2385. 

2385.  2385.  2382 

.  2374. 

2368. 

2369. 

2378. 

2383. 

1 

2377. 

2368. 

2367. 

2378 

2371.  2375. 

2376. 

2369. 

2363. 

2365. 

2365 

25 

22 

2848 

2377. 

2368  .  2367  .  2378 

.  2371. 

2375. 

2376. 

2369. 

2363. 

1 

2367. 

2366. 

2369. 

2372 

2373.  2373. 

2369. 

2378. 

2374. 

2376. 

25 

23 

2848 

2367. 

2366.  2369.  2372 

.  2373. 

2373. 

2369. 

2378. 

2374. 

2376 

1  2377. 

2381. 

2387, 

8391 

2388. 

2382. 

8379. 

2376. 

2374. 

2371. 

as 

NO.  OF  CAL  BUS 
15 

A 

24 

25 

2048 

0 

2377. 
NO.  OF 

2381.  2387.  2391.  2388.  2382 

DATA  BUS  25 

.  2379. 

2376. 

2374. 

2371. 

161 

2048 

2168. 

2165. 

2162* 

2163. 

2162. 

8160. 

2158. 

2158. 

2163. 

2167. 

HI 

2048 

3184. 

3188. 

3196. 

3191. 

3194. 

3196. 

3192. 

3186. 

3181. 

3181. 

Hi 

16 

A 

2048 

85 

2377. 

2381. 

2387. 

2391. 

2388. 

2382. 

2379. 

2376. 

2374. 

2371. 

201 

2048 

2042. 

2037. 

2635. 

2632. 

2629. 

2028. 

2027. 

2036. 

2044. 

2045. 

201 

2048 

3029. 

3022. 

3616. 

3618. 

3623. 

3026. 

3023. 

3025. 

3031. 

3036. 

201 

17 

A 

2048 

25 

2031. 

2036. 

2646. 

2644. 

2646. 

2040. 

2042. 

2042. 

2045. 

2041. 

301 

2048 

2726. 

2725. 

2719. 

2722. 

2724. 

2722. 

2720. 

2727. 

2732. 

2735. 

301 

2048 

1204. 

1197. 

1261. 

1216. 

1213. 

1217. 

1218. 

1216. 

1206. 

1199. 

301 

18 

A 

2048 

25 

3412. 

3413. 

3463. 

3398. 

3464. 

3411. 

3410. 

3403. 

3400. 

3411. 

4> 

401 

2048 

2565. 

2569. 

2575. 

2576. 

2569. 

2559. 

2552. 

2552. 

2559. 

2564. 

O' 

401 

2048 

1183. 

1181. 

1186. 

1178. 

1179. 

1188. 

1195. 

1195. 

1188. 

1177. 

401 

2048 

3405. 

3402. 

3465. 

3416. 

3416. 

3401. 

3399. 

3403. 

3405. 

3404. 

T  yic 

/framf;--  25.«C0 

/TIC-  6.250 


4000 


APPENDIX  D 

IMPLEMENTATION  OF  STEP  3 
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APPENDIX  D.l.a 

Listing  of  MFA  File  ADTAPE3/UN=B00TS 
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APPENDIX  D.l.a 


Listing  of  MFA  File  ADTAPE3/UN=B00TS 


a.  Job  Control  Language 

BOOTS, STMFZ,P6,T30. 

ACCOUNT, PDxxx. 

BEGIN , ATTACH , PLOTLIB . 

ATTACH , TAPE 1 , LOLAD , ID=BOOTS . 
REQUEST, TAPE3,*PF. 

REQUE ST , TAPE 1 3 , * PF . 

FTN,R=0. 

MAP, OFF. 

LGO. 

CATALOC , TAPE 3 , LOLAED , ID=BOOTS . 
CATALOG , TAPE 1 3 , PLT , ID=BOOTS . 
BEGIN , PLOT , CALCOMP , TAPE 1 3 . 
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APPENDIX  D.l.b 
Program 
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PROMMI  EDITADC  |HWr.0UTWT,T<«5-IHPUT,T<#M-0UfTPUT.TMCl.T^f3. 
*C0m0N9TEn*  <  8S«* ).  COUNT  ( 83M«  )#  MTA  ( 83«M  ).ZPLt  ICMT,  NTS  ( 84 ) , 

* *5E«i35!ncma  .ImiTWiEfTdtInPT.TXNC,  if*aco5,  isaAu ) 


IDIH-83800 
UNITE <8* 1800) 

68  CALL  ND&TAC1) 

IPL-1 

DO  500  I-l.NCH 

ICNT-I 

NPT-0 

DO  100  N»1#NCHA 
CALL  RDDTAC8) 

CALL  EDITCN) 

188  CONTINUE  _  _ 

IF<NPT.E0.8)C0  TO  580 
CALL  URTTPE 
588  CONTINUE 
CO  TO  58 
1888  FORMAT (1H1 > 

END 


1 


SUBROUTINE  EDITCN) 

CONNON  TEMP(8588).C0UNTC83888),DATA(83888)»IPl#ICNT»NT8<84) 


i 


6 

1# 

15 

2* 

25 

3t 

35 


commotT  TwfisiS )  {coSfcgjM ) ,  0J^» ; ?2J^  >  fci^.icrrr,  MTS  ( 84  > , 

t  NCH,NU,IRD,NCHA.lDIN,NAHEF,DT,NPT,TINE,IFRA<S>,ISAH<8> 
IF(ICODE.NE.l)  Q&  TO  8M 
READtl >NAMEF  _  _ 

IFCEOF(l).NE.O. )  STOP 
REM  ( 1 )  IRD  *  DT.NCH.NCMA 
URITEC6. 1  >NAMEF,HCH,DT,NCHA 
IF<IRD.LT.1#00A>URITE<S,4) 

IFCIRD.QE.1MA0>URITE<6,5) 

READ  <  S  ,  3 )  IFRA  ,  ISAH.IHSG 
IF<E0FC5).NE..0>ST0P 
DO  1M  1*1.24 
NTS(I >*2H 
1M  CONTINUE 

IFdNSC.NE.#>READ<S,6XNTS<I>, 1*13,24) 

TIHE-.S 

IF  < ISAMd ).EQ.S  >ISAM( 1  )*1 
IF(ISAH(2).EQ.0)ISAH(8)*1 
URITEC6,3>IFRA,ISAH,INSC 
DT*DT*FL0AT(ISAM(2>) 

RETURN 

200  IFCIC0DE.NE.2)  GO  TO  30# 

READ(l)  HU, <TEHPCI>,I*l»NU) 

IFCEOFC 1 ).NE.O. )  STOP 
RETURN 

3M  URITE(6,2  >  I  CODE 

1  FoB5t<//.'  FILE  NAHE  '.AlO.lSX.'NO.  OF  CHANNELS, •  l,JS/_ 


r  0  r  ibc.  iwwifc  f  * 

t  '  TINE  INC RENE NT, SEC  ■ 

*  15// > 

2  FORNATC '  RDDATA  ERROR,  ICODE  ■  Mil) 

3  FORMAT (915) 

4  FORMAT ( *  DATA  CHANNEL'  ) 

5  FORMAT ( '  CALIBRATION  CHANNEL' ) 

6  FORMAT ( 12A2 ) 

END 


.E12.S,'  NO.  OF  DATA  SUBSETS 


58 


1 

s 

It 

IS 

29 

2S 

39 

35 

49 

45 

59 

55 


CoSSSTTCHPCaSi*  Jc8WS3999 >,DATA<23999>, IPL, ICNT,NTS<24  >, 
tJSuJu!  IRD,HCHrt,  J0W.N9«CF,»T,r#*T,  TIME,  159^(2 ) 

DINENSlfo  LAiEL<4)JTlTLE(3) 

IF  < IPLOT . OT . I )  00  TO  199 

IPLOT-2 

XB-1.7S 

VB-.4 

XS-699. 

DX-S99. 

XNI-.9 

XAX-NPT/599*1 

XNA-XSXXAX 

VAX-4. 

XPA6E-XAX+2.2 

VPAQE-4.75 

FACT-1. 

I UNIT- 13 

LABEL! 1 >-19H  BOOTS  399 
LABEl(2)-19N  6121 
LABEL(3)-19HE0ITTE0  DA 
LABEL(4)-19HTA  -  L0LA2 
NODE-1 

CALL  PLTBEfi(XPAGE , YPAGE, FACT, IUNIT, LABEL ) 

YNI-9. 

VNA-5999. 

DV-1999. 

VS-(VNA-vm>/VAX 

199  CALL  PLTSCA<XB.V1,XNI,VNI,XS,YS) 

CALL  PLTyND(XNI,XRA,VNI,VNA) 

CALL  PLTAX5( DX, DV , XNI , XNA, YNI , VNA, 4 ) 

CALL  LADELA<DX,DY.XNI.XNA,YNI,YNA,1.,1.> 

ENC0DEC39, 1999, ITITLE)  IRD1 
1999  FORMAT (3H  ID,I19,8H  PL0T,8X,1H>> 

TX-XNI 

TV-VNA+.9SBVS 


CHT-.i 

CALL  PLTSVH(CHT,ITITLE(1),9.,TX,TV) 
TIC-DXSDTX1999. 

TIN-DT*FLOAT(NPT  >*1999. 
TV-YNA+.29*YS 

ENCODE ( 12,2999, ITITLE )NANEF 
2999  FORMAT ( 1H  ,A19,1H>  ) 

CALL  PLTSYM(CHT, ITITLE, 9., TX,TY1 
TX-XRI-XSX1.5S 


TY-YNA-VSX.3 

ITITLE ( 1 )-19MTIHE,  NS> 

CALL  PLTSYN<CHT, ITITLE, 9., TX,TY> 

TV- VNA- VS*. 45 
ENCODE ( 16, 3999 .ITITLE )T1H 
3999  FORMAT (BN  /FRAHE-,F7.3,1H> > 

CALL  PLTSVNICHT, ITITLE,9.,TX,TV> 
TY-VNA-VSX.6 

ENCODE! 16,4999, ITITLE >TIC 
4999  FORMAT! 8N  /TIC-  ,F7.3,1H>> 

CALL  PLT$YN<CHT.ITITL£,9.,TX.TV> 

CALL  PLTDTS(N0DE,9,C0UNT(1>,MTA(1),NPT,9) 

CALL  PLTPOE 

RETURN 

END 
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COM^TENP(SS^.  COUHT  ( 83tM ).  MT«(  IMM ).  IPL.  ICHT.  NT$(  84 ) , 
*  NcJmJuIiRdI  NCHA.  blN,)*NEF,DT.NPT,TINE,  IF*A<0>,  ISAACS ) 

SO  5  I-1,XDIN 
COUNT (I >■! 

6  DEC0Dl!¥l0,8000,NANEF>(NTS(I),I-l,S) 

ENCODE (  6, 3000, NOUN  )IAD 

DECODE (  6,8000,NDUN  HNTSCI  ),  I-t#  It) 

ENCODE <8,4000,NTS< 11 >  >ICNT 
ENCODE (  $,4000,NTS<12)>NCH 
URXTE(3)NTS 
URITEO)  NPT 

URITEO  WINE, DT  _ 

URITEO, 1000 (NTS, NPT, TINE, DT 
CALL  PLAT  ( IPL, IRD) 

URITEO)  (DATAtJ ), J»1,NPT) _ 

URITEO,SOOO  )<DATA(  J >, J*l,80) 

RETURN 

1000  FORNAT < '  NTS  -  ;,84A8,'  NPT  -  ',110,'  START  TINE  •  ', 

t  E18.S,'  DT  -  ',E18.S) 

2000  FORNAT (SA2) 

3000  FORNAT ( IS ) 

4000  FORNAT (18) 

SOOO  FORNATOOFIO.O) 

END 
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APPENDIX  D.I. C 
Input  Card  Images 
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3 

15 

16 

16 

16 

16 

58 

1 

0 

4 

5 

5 

6 

6 

27 

56 

2 

0 

2 

13 

14 

14 

14 

18 

58 

1 

0 

5 

6 

6 

7 

7 

28 

50 

2 

6 

3 

4 

4 

5 

5 

26 

56 

2 

0 

4 

15 

17 

17 

17 

19 

56 

1 

9 

63 


APPENDIX  D.2 
Sample  Output 
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FILE  NAME  L2R47 
TIME  INCREMENT, SEC 


4 


25 


NO.  Of  CHANNELS  • 
.12S66E-64  NO.  Of  DATA  SUBSETS  • 


DATA  CHANNEL 

3  15  16  16 

NTS  •  L2R47 

2121.  2123. 

2142.  2146. 

NTS  •  L2R47 

2036.  2636. 

8036.  2043. 

NTS  •  L2R47 

8624.  2613. 

8679.  2696. 

NTS  ■  L2R47 

8433.  2443. 

8499.  2562. 


16  18  56  1 

1  1  4 

8186.  2136. 

2141.  2152. 

12  4 

8638.  2639. 

8636.  2636. 

1  3  4 

2626.  2688. 

2678  .  2669. 

1  4  4 

2427.  2446. 

2S67.  2S17. 


6 

NPT  • 

2134.  2136. 

8156.  2146. 

NPT  - 

8647.  2639. 

2636.  2641. 

HPT  • 

2651.  2632. 

8766.  2693. 

NPT  • 

2445.  2453. 

2518.  2S35. 


O' 


HO.  Of  CHANNELS  -  1 

.  12566E-64  HO.  Of  DATA  SUBSETS  •  28 


FILE  NANE  DP47 
TINE  INCREMENT, SEC  - 


DATA  CHANNEL 

4  5  5  6 

NTS  •  DP47 

2636.  2629. 

2636.  8645. 


6  27  56  2 

2  1  1 

2646.  2636. 

2633.  2645. 


6 

NPT  ■ 

2629.  2638. 

2638.  2636. 


FILE  NAME  L8R46 

NO.  OF 

CHANNELS  ■ 

4 

TIME  INCREMENT, SEC  - 

. 12566E- 

■64  NO.  OF 

DATA  SUBSETS 

■ 

DATA  CHANNEL 

2  13 

NTS  >  L2R46 

14  14 

14  18 

3  1  4 

56  1 

0 

NPT 

2611. 

2607. 

8019. 

2616. 

2004. 

2017 

2664. 

NTS  -  L2R46 

2011. 

2004. 

3  2  4 

2661. 

2006. 

2016 

NPT 

2356. 

2336. 

8328. 

2336. 

2338. 

2358 

2346. 

NTS  •  L2R46 

2350. 

2342. 

3  3  4 

2352. 

2336. 

2341 

NPT 

2965. 

2978. 

2978. 

2986. 

2973. 

2991 

2974. 

NTS  •  L2R46 

2971. 

2957. 

3  4  4 

2969. 

2976. 

2976 

NPT 

2845. 

2833. 

2856. 

2849. 

2830. 

2850 

2839. 

2817. 

2843. 

2817. 

2833. 

2837 

4588  START  TIME  -  6. 
2136.  2139.  2133. 
2151.  2152.  2149. 

4588  START  TIME  -  6. 
2633.  8636.  8642. 
2633.  2643.  2643. 

4588  START  TINE  -  6. 
2644.  2646.  2651. 
2686.  2689.  2676. 

4588  START  TIME  -  6. 
2456.  2462.  2482. 
2513.  2S27.  2529. 


22569  START  TIME  -  6. 
2644.  2632.  2647. 
2641.  2635.  2641. 


8643  START  TIME  -  6. 
2661.  2615.  2669. 
2613.  2669.  2616. 

8643  START  TIME  -  6. 
2346.  2351.  2352. 
2337.  2344.  2329. 

8643  START  TIME  -  6. 
2972.  2967.  2972. 
2989.  2984.  2973. 

8643  START  TIME  -  6. 
2815.  2839.  2835. 
2832.  2864.  2847. 


DT  • 

.125666-64 

2133. 

2157. 

DT  • 

. 125666-64 

2033. 

2041. 

DT  ■ 

. 12566E-64 

2653. 

2691. 

DT  • 

. 12566E-64 

2496. 

2531. 

DT  ■  .256666-64 

2633. 

2633. 


DT  ■ 

. 125666-64 

2666. 

2629. 

DT  • 

. 12S66E-64 

2356. 

2318. 

DT  • 

. 125666-64 

2966. 

2985. 

DT  • 

. 12566E-64 

2817. 

2851. 

CviO 


APPENDIX  E 

IMPLEMENTATION  OF  STEP  4 


69 


APPENDIX  E.l 

Listing  of  MFA  File  ADTAPE4 /UN=BOOTS 


71 


APPENDIX  E.l 


Listing  of  MFA  File  ADTAPE4/UN=B00TS 


a.  Job  Control  Language 


BOOTS , STMFZ , P6 ,MS300000 . 
ACCOUNT , PDxxx . 

REQUEST , NEWPL , *PF . 

REQUEST, BIN, *PF. 

BEGIN , ATTACH , PLOTL IB . 

ATTACH , OLDPL , ADENGR, ID=BOOTS . 
UPDATE, F,N. 

FTN , I=COMP ILE , L=0 . 

ATTACH, TAPE 1 ,LOLAED, ID=BOOTS. 
LGO. 

EXIT, U . 

BEGIN , PLOT , CALCOMP , TAPE 1 3 . 
*EOR 


b.  UPDATE  Corrections 

* I DENT  HP ADI 
*EOR 


c .  Input  Card  Images 

14  1.  1000  .0  0  1 

$CONVAR  NX=50,NY( 1)=500,NY(2)=70,NY(3)=150,NY(4)=230,NY(5)=320, 
NY( 6)=400 , ITZ=1000 , IBSE=3700$ 


73 


APPENDIX  E.2 

Listing  of  MFZ  File  ADENGR, ID=BOOTS 


75 


1 


*^J 

*^J 


5 

1# 

15 

25 


25 

30 

35 


45 

55 

55 


C 

C 

c 

c 

c 

c 

c 


c 

c 

c 

c 


PROGRAM  DATA  ( INPUT, OUTPUT, TEAR, TAPE  1  ,TAPE2-TEMP, TAPES- INPUT, 
rr  APE6-0UTPUT ,  TAPE  1 3 ,  TAPE3 ,  TAPE4 ) 

FOR  THIS  PROGRAM  BROAD  HAS  KEN  MODIFIED  SO  THAT 
THE  ENTIRE  DATA  SET  IS  READ  IN  AT  ONE  TIME* 

N  CANNOT  EXCEED  23000; 

SINGLE  INTEGRATION  CAN  K  DONE  ON  THE  ENTIRE  DATA  SET; 

DOUBLE  INTEGRATION  CAN  BE  DONE  IF  N  DOES  NOT  EXCEED  l2SOO. 

COMMON  RC2S000  >,TC23000>,PC23008>, NOS, LP,LABEL<4),NYC6),Y<6>,  IOCS) 

1  ,SBL,I,TZ,JP,IX,DELT,IDZ,NOP<10>,LLCtB,SA,SB,DTIHE,NJ<10KTSTART, 

2  SC,  SEP.  XAX,  XFAC,  YFAC,  TITC30),  ISU,  KP ,  VI,  XB,  NTS(24>,ITZ, 

3  JPICSO ), iZ, ITM, NOSR,N, IPN0C5  ),SS,IC$TS1  ,ICSTS2,  VAX,  JSU,KU 
DIMENSION  PIC1),PI2(1> 

EQUIVALENCE  (RC 1 >,PI<1 )  >, (R( 12501 >,PI2C 1  ) ) 

DIMENSION  SO),  RU<6>,  AF(6),  SIGC3),  TGC6),  X(6),  E0C6, 

1  3).  FO,4),  SXC6) 

HAMEL  IST/CONVAR/SKP ,  XAX ,  VAX  ,  XF  AC ,  YFAC  ,  I TM,  NX ,  NY ,  MOP  ,  SS  ,  I TZ ,  NS  ,  I  BSE 

LABELC1 >-10H  BOOTS 

LAKLC2>-10H390  6121 

LABELC3)-10H  A/D  DATA 

LABEL ( 4 ) - ItHREDUCTION 

UR  I  TEC  6,  MO) 

TITC26  )■  1K> 

ITM-23800 
ITZ-  950 
LLCS-  1 
MOPC 1  )■  1 
N0PC2)-  3 
M0PC5  )•  • 

NX-40 
NVCD-16 
NY (2 )*100 
NYO  )-240 
NVC4  )-350 
NY(S)-4M 
NYC  6) -600 
SJCP-159* 

SS-  1. 

XAX-7.5 
XFAC-  I, 

YAX-  5* 

YFAC-  1. 

NS-6 
I BSE -500 

READING  DATA 

XB-2. 

VB-1.5 

ICSTSl-0 

KSTS3-0 

20  IFCNUC«GT#0  >G0  TO  21 

READ(S,7O0>IPN0<2),IPN0C3>,B,IXS,DLTN,IDZS,NLIST,IFID,NU 
IF  CEOF(S)  .HE.  0.)  STOP 
IF  CNLIST  .NE.  0)  READC 5, CONUAR ) 


MAIN 

C19 

MAIN 

HPAD 

HPAD 

HPAD 

HPAD 

HPAD 

HPAD 

HPAD4 

HPAD 

COM 

HPAD 

HPAD 

HPAD 

MAIN 

MAIN 

HPAD4 

C2 

HPAD3 

HPAD 

C3 

Cl 

MAIN 

HPAD3 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

MAIN 

C3 

MAIN 

MAIN 

MAIN 

CAL2 

CAL  4 

MAIN 

MAIN 

MAIN 

MAIN 

C3 

C3 

MAIN 

MAIN 

HPAD2 

HPAD2 

MAIN 

MAIN 


8 

2 

4 

7 

8 
9 

10 

11 

12 

1 

3 

4 

4 

5 

6 
€ 
7 
2 
3 
1 

13 

2 

1 

9 

2 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

3 

26 

27 

28 
2 
2 

29 

30 

31 

32 

4 

5 

35 

36 
1 
2 

38 

39 


66 


C 

C 

c 


oo 


65 

76 

75 

86 

85 

96 


95 


166 


165 


116 


21  NU-NUC-1 

HPAD2 

3 

DO  25  ICT-1,16 

C19 

3 

25  NJOCT  )-6 

C19 

4 

MIN 

46 

INTERPRETING  OPTIONS  FROR  3-DIGIT  CODE 

MIN 

41 

MIN 

42 

IF(IPN0(2).EQ.6*AND.IPN0<3).EQ.6)G0  TO  666 

HPAD 

15 

11*2 

HPAD 

16 

IF  ( IPNO( II ) *NE* 1 )  GO  TO  126 

MIN 

81 

IDZ-1 

MIN 

82 

IX-566 

MIN 

83 

GO  TO  266 

MIN 

84 

126  IF  ( IPNOni  >«NE»2 )  GO  TO  136 

MIN 

85 

IDZ-1 

MIN 

86 

IX-2666 

MIN 

87 

GO  TO  266 

MIN 

88 

136  IF  ( IPNOCII )  .NE.3  )  GO  TO  146 

RAIN 

89 

IDZ-1 

MIN 

96 

IX-3666 

MIN 

91 

GO  TO  266 

MIN 

92 

146  IF  (IPNOCII). NE. 4)  GO  TO  156 

MIN 

93 

IDZ-3 

MIN 

94 

IX-566 

MIN 

95 

GO  TO  266 

MIN 

96 

156  IF  (IPNOCII ).NE.5)  GO  TO  166 

MIN 

97 

IDZ-3 

MIN 

98 

IX-2666 

MIN 

99 

GO  TO  266 

MIN 

166 

166  IF  ( IPNOCII ).NE.6 )  GO  TO  176 

MIN 

161 

IDZ-3 

MIN 

162 

IX-3666 

MIN 

163 

GO  TO  266 

MIN 

164 

176  IF  ( IPNO( II ) «NE« 7 )  GO  TO  186 

MIN 

165 

IDZ-5 

MIN 

166 

IX-566 

MIN 

167 

GO  TO  266 

MIN 

168 

186  IF  ( IPNOCII )«NE«8 )  GO  TO  196 

MIN 

169 

IDZ-5 

MIN 

116 

IX-2666 

MIN 

111 

CO  TO  266 

MIN 

112 

196  IF  (IPNOC II >.NE.9 )  GO  TO  676 

MIN 

113 

IDZ-5 

MIN 

114 

IX-3666 

MIN 

115 

266  II-3 

MIN 

116 

IF  ( IPNOCII ).NE. 1 )  GO  TO  216 

MIN 

117 

KP-6 

MIN 

118 

N0P<3)-6 

RAIN 

119 

N0P<4)-6 

RAIN 

00  TO  296 

RAIN 

216  IF  ( IPNOCII ) «NE»2 )  GO  TO  226 

RAIN 

|HT7j 

KP-6 

MIN 

123 

ROPOJ-2 

MIN 

124 

R0PC4  )-6 

MIN 

125 

GO  TO  296 

MIN 

126 

226  IF  (IPN0< II ) *NE*3 )  GO  TO  236 

MIN 

127 

ICP-6 

MIN 

■iij 

N0P<3)-2 

MIN 

129 

115 


185 


185 


139 


135 


149 


145 


159 


155 


169 


165 


C 

C 

C 


N0P<4>-8 
60  TO  899 

839  IF  ( IPNO( II ) *NE.4  )  60  TO  849 
1CP"3 
N0P<3>-9 
N0P<4)-9 
60  TO  899 

849  IF  (IPNO(II  )»N£*5 )  60  TO  859 
ICP-8 
N0P(3>-8 
N0P<4>-9 
60  TO  899 

859  IF  ( IPHO( II )«NE*6 )  60  TO  869 
KP-8 
N0P(3)*2 
N0PC4)*8 
60  TO  899 

869  IF  <IPN0<II)*NE.7>  60  TO  879 
KP-1 

H0P(3  )*9 
HOP  (4  )*9 
60  TO  899 

879  IF  (IPNOCII )«NE»8 )  60  TO  889 
KP*1 

N0P<3)-8 
N0P<4>-9 
60  TO  899 

889  IF  ( IPMO< II ) *NE»9  )  60  TO  679 
KP-1 

N0P<3>*2 
N0P(4>»2 
60  TO  899 

899  IF  (IXS.NE.9)  IX-IXS 
IF  (IDZ$'NE«9>  IDZ-IDZS 
KU*9 
isu-9 
JSU-9 

READING  FRON  TAPE 

NOS-1 
NOSR-9 
N-l 
IZ-9 

DO  649  IP- 1, NOS 
LLC-LLCS 

IF  (ICU»E6«9 )  I»3 
CALL  DATAIN 

TZ*DELTtl999  •  SFLOAT  (IX  J-DLTIHTSTART 
DTINE-FLOAT  ( IDZ  XDELTS1999* 

J*1 

IF  (ICU.E6.19)  60  TO  549 
IF  (LP.GT.l)  60  TO  569 
URITE(6,719HNTS(K),K-1,12> 


179 


MAIN 

13* 

MAIN 

131 

MAIN 

138 

MAIN 

133 

MAIN 

134 

MAIN 

135 

MAIN 

136 

MAIN 

13? 

MAIN 

138 

NAIN 

139 

MAIN 

149 

NAIN 

141 

NAIN 

142 

NAIN 

143 

NAIN 

144 

NAIN 

145 

NAIN 

146 

NAIN 

147 

NAIN 

148 

NAIN 

149 

NAIN 

156 

NAIN 

151 

NAIN 

152 

NAIN 

153 

NAIN 

154 

NAIN 

155 

NAIN 

156 

NAIN 

157 

NAIN 

158 

NAIN 

159 

NAIN 

16# 

NAIN 

161 

NAIN 

162 

NAIN 

164 

NAIN 

165 

NAIN 

16? 

NAIN 

168 

NAIN 

169 

NAIN 

17# 

NAIN 

171 

HP  AD 

17 

HPAD 

18 

HPAD 

19 

C19 

5 

NAIN 

178 

NAIN 

179 

NAIN 

18# 

HPAD 

21 

HPAD 

22 

HPAD 

23 

NAIN 

194 

NAIN 

195 

NAIN 

196 

HPAD 

24 

NAIN 

198 

NAIN 

199 

NAIN 

2## 

C 

C 

C 


PREDICT  METHOD  OF  CALCULATING  CAL  I  DRAT  I  ON  STEPS 


00 

o 


SUN*#* 

2  *NV ( J  ) 

IPP*1 

IOC-1 

SUNN*#* 

IT(A#S<R(  I  )-R(  1-1 ) )  .LT.SICP )  GO  TO  4## 

IF  Aftl  R<I>-R<1-2)>*LT*$ICP)  00  TO  39# 

IF  <A#S<R<I*1>-«<I-1)>.GT.SKP)  CO  TO  37# 

IF  ™»I<R<I+2>-mi-l>>.CT.SICP>  CO  TO  380 

i»i*a 

CO  TO  400 

IF  (ICK.GE.IO)  CO  TO  650 

ICK-ICK+l 

I-I+3 

CO  TO  360 

I-I+l 

SUH*SUN+R(I> 

suhh-suhn+r<I)s*2 

IPP-IPP+1 

1-1*1 

IF  (IPP.EQ.NX+1)  GO  TO  410 
GO  TO  360 

X(J)-SUN/FLOAT(NX) 

SX(  J ) -SORT( SUNN/FLOAT ( NX )- ( X ( J  )»2 ) ) 
IF(J.EQ.NS)  GO  TO  420 
J-J+l 
GO  TO  350 

FORMING  EQUATIONS  AND  GENLSQ 


V(1 >-0. 

V(2)-1.SSS 

V(3>-2.SSS 

V(4)-3.X$S 

V(5)-4.S$S 

V(6)-5.SSS 

SBL-X(l) 

WRITE(6»720 )  (X(n),A-l,NS). (SX(N).N-l.NS) 

DO  430  N-1*NS 
X(N)-X(N)-SBL 
DO  440  L*l*6 
EQ(L*1)-1. 

EQ(L,2)-X(L) 

EQ(L*3 )-X(L )St2 
WRITE  (6*730) 

DO  450  L-1,NS 

WRITE  (6*740)  EQU,  l  ).EQ(L,2)*EQ(L*3 ),V(L) 

IF  (A1S(X(1)-X(2)).IE.20. )  CO  TO  460  , 

CALL  GENLSQ  (EQ.6,V,NS,F*3,3,S.RU,AF.ERNS,SIC,T6,DeT,0 ) 
SA-S(t) 

SB-S(2) 

SC-S(3) 

GO  TO  470 
SA*0. 

Sl-i. 

sc-o. 


225 


B-l. 

NAIN 

259 

230 

SIG(l)-0.0 

RAIN 

200 

$IG(2)-0.0 

NAIN 

201 

SIG(3)-0.0 

NAIN 

202 

470  URITE  (6,750) 

NAIN 

263 

URITE (6, 760 )NTS( 10),NT$(11  ),NTS(  12 ),B,SA,$B,$C,SIG(1  ),SIG(2  ),  HPAD 

26 

235 

*  SIGO) 

HP  AD 

26 

0 

NAIN 

265 

C  SEARCHING  FOR  FIDUCIAL 

NAIN 

266 

C 

NAIN 

267 

i-i+itz 

NAIN 

268 

24# 

IF(IFID*NE.0)G0  TO  525 

HPAD 

27 

480  IF  ( AiS(R(  I  )-R(  1-1 )  ).LT.SKP )  GO  TO  500 

NAIN 

269 

IF  (AB$(R(I)-R(I-2)).LT.SKP)  GO  TO  500 

NAIN 

270 

IF  (ABS(R( 1*1 )-R(I-l )  KGT.SKP)  GO  TO  490 

NAIN 

271 

GO  TO  510 

NAIN 

272 

245 

490  IF  (ABS(R(I+2)-RCI-l)).GT.$JCP)  GO  TO  520 

NAIN 

273 

I-I+3 

NAIN 

274 

GO  TO  480 

NAIN 

275 

500  I-I+l 

NAIN 

276 

GO  TO  480 

NAIN 

277 

250 

510  1*1+2 

NAIN 

278 

GO  TO  480 

NAIN 

279 

520  URITE  (6,770)  1 

NAIN 

280 

525  SUN-,0 

HPAD 

S3 

IFS-I+IBSE 

CAL4 

255 

DO  530  KF-1,100 

NAIN 

283 

$UN-SUH+R(IFS) 

NAIN 

284 

IFS-IFS+l 

NAIN 

285 

530  CONTINUE 

NAIN 

286 

SOL -SUN/ 100* 

NAIN 

287 

260 

URITE  (6,780)  SBL 

NAIN 

288 

I-I+IX 

NAIN 

289 

540  IF( I ,LT .  JP )G0  TO  560 

HPAD1 

1 

JPI(LP)-0 

NAIN 

291 

KU-10 

NAIN 

292 

265 

GO  TO  680 

HPAD1 

2 

C 

NAIN 

296 

C  CALLING  SUBROUTINES 

NAIN 

297 

C 

NAIN 

298 

560  NIP-NOP(LLC) 

NAIN 

299 

270 

IF  (NIP, EG. 0)  GO  TO  630 

NAIN 

300 

GO  TO  (580,590,570),  NIP 

NAIN 

301 

570  IF  (KP.EQ.O)  GO  TO  620 

NAIN 

302 

GO  TO  (600,610),  KP 

NAIN 

303 

580  CALL  PRESS 

NAIN 

304 

275 

LLC-LLC+1 

NAIN 

305 

GO  TO  560 

NAIN 

300 

590  CALL  INTO 

NAIN 

307 

LLC-LLC+1 

NAIN 

308 

GO  TO  560 

NAIN 

309 

280 

600  URITE(6,601) 

HPAD 

29 

601  FORNAT(  *  THE  DECK  PL0T1  HAS  BEEN  YANKED' ) 

HPAD 

30 

LLC-LLC+1 

NAIN 

311 

GO  TO  560 

NAIN 

312 

610  ISU-ISU+1 

NAIN 

313 

285 

CALL  PL0T2 

NAIN 

314 

299 


296 


3N 


395 


319 


315 


329 


325 


C 

C 

c 


629  UC*UC+1 
60  TO  569 
639  ICU-9 
649  CONTINUE 

CALL  DATASOCN) 

WRITE (6, 799  )NTS( 19 ),NTS( 11  ),NTS(  12),KN 
WRITE  (6,899 ) 

REWIND  2 
REWIND  4 
GO  TO  29 

ERROR  PRINTS 

659  WRITE  (6,819)  R(I >,R(I+1  >,R(I*2) 

STOP 

669  READ  (5,829)  IQ 
CALL  SKIP 
GO  TO  29 

679  WRITE  (6,839)  II,IPN0(I1) 

STOP 

689  WRITE (6, 849 )I, JP 
GO  TO  29 

799  F0RMAT<3X,2I1,F19#3, 1 19, F 19 •9,915 ) 

719  F0RMAT(///#4H  ID-,12A2) 

729  FORHAT  (6H  X'S  ,6(F19.3,2X)  /  SX-  6(F19.3,  2X)) 

739  FORMAT  ( 1H9,26X,  18H0RIGINAL  EQUATIONS) 

749  FORMAT  <1H  ,4<E1S*8,2X ) ) 

759  FORMAT  ( 1H9,27HCHANNEL  GAGE  CAL  IB  CONSTANT,  6X,  2HSA,  1 IX,  2HSB,  1 IX,  2H 
1$C,9X,6HSIG(  1  ),7X,6HSIG(2 ),7X,6H$IG<3  > ) 

769  F0RMAT(2X,3A2,3X,F12.4,5X,6(E12.5,1X)  ) 

779  FORMAT  <  1H9, 19HINDEX  TO  FIDUCIAL  -.119) 

789  FORMAT  ( '9BASELINE  AVERAGE  -  ',F9.2) 

799  FORMAT <  1H  ,27HTHIS  IS  THE  END  Of  CHANNEL  ,3A2,'  - 
t'  DATA  POINTS  SAVED' ) 

899  FORMAT  ( 1H1 ) 

819  FORMAT  (1H  ,36HTHERE  IS  A  LOT  Of  SOMETHING  GOING  0N,3X,3(F19-3,3X ) 
1) 

829  FORMAT  (15) 

839  FORMAT  (31H9THIS  OPTION  IS  INCORRECT,  II  -,IS,3X,3X,8H0PTI0H  -,I5> 
849  FORMAT( '  I.GE.JP,  I  •  ',119, '  JP  •  ',119) 

END 


,119, 


MAIN 

315 

MAIN 

318 

MAIN 

317 

MAIN 

318 

HPAD 

31 

HPAD 

38 

MAIN 

388 

HPAD 

33 

C21 

1 

MAIN 

388 

MAIN 

383 

MAIN 

324 

MAIN 

386 

MAIN 

386 

MAIN 

327 

MAIN 

328 

HPAD4 

3 

HAIM 

338 

MAIN 

331 

MAIN 

vmm 

HPADl 

3 

HPAD1 

4 

MAIN 

333 

HPAD 

35 

HPAD4 

4 

MAIN 

338 

MAIN 

337 

MAIN 

338 

MAIN 

339 

MAIN 

348 

HPAD 

37 

MAIN 

342 

MAIN 

343 

HPAD 

38 

HPAD 

39 

MAIN 

345 

MAIN 

348 

MAIN 

347 

MAIN 

348 

MAIN 

349 

HPAD1 

5 

MAIN 

358 

1 

SUBROUTINE  DATAF  _  . 

DAT  AF 

a 

COMMON  RC2S990>,T<23000>,P<23990),N0S,LP,tAIEt<4>,NV<6>.V<6).IQ<S> 

1  ,Slt,I,TZ,JP,IX,DElT,IDZ.H0Pd9).llC,B.SA,SI.DTIHE.NJd9>,TSTART. 

2  SC,  SKP,  XAX,  XFAC,  VFAC,  TITOi),  ISU,  KP,  VI.  XI.  NTS<24).ITZ. 

HPAD4 

i 

HPAD 

3 

CON 

4 

5 

3  JPH59),  12, ITM,NOSR,N, IPN0(5 ),SS,KSTS1 ,KSTS2,VAX, JSU,KU 

HPAD 

4 

DIMENSION  PI(1 ),PI2( 1 ) 

HPAD 

5 

EOUIUALENCE  (Rd  ),PI(1)),(R( 12591 ),PI2( 1 ) ) 

HPAD 

6 

IF( JSU.EQ.l )C0  TO  29 

C19 

7 

URITE<4)NTS 

C19 

8 

19 

URITE(4)T(1 ),DTINE 

C19 

9 

URITE(4)N 

C19 

19 

URITE(4)CP(IT).IT-1.N) 

C19 

11 

NJ(LP)-N 

cat 

1 

60  TO  39 

C19 

13 

15 

29  NJ(LP>-9 

C19 

14 

39  RETURN 

C19 

15 

00 

u> 

END 

DATAF 

11 

1 

C 

HPAD 

49 

SUBROUTINE  DATAIN 

HPAD 

41 

C 

HPAD 

48 

COMMON  R (35999 ), T< 23999 ),P(23999), NOS, LP,LAIEL(4),NV(6), V(6),IQ(5) 

1  ,SBL,I,T2, JP, IX.DELT, IDZ,NOP( 19),LLC,B,SA,SI,DTINE,NJ(  19),T$TART, 

2  SC.  SKP,  XAX,  XFAC,  VFAC,  TIT(39).  ISU.  KP.  VI,  XI,  NTS<34).ITZ, 

HPAD4 

1 

5 

HPAD 

3 

CON 

4 

3  JPK59 ).IZ, ITN,NOSR,N,IPNO(S ),SS.KSTS1,KSTS2,YAX, JSU.KU 

HPAD 

4 

DIMENSION  PI(1),PI2(1) 

HPAD 

5 

EOUIUALENCE  (R( 1  ),PI( 1 ) ), <R( 12S91 ),PI2( 1 ) ) 

HPAD 

€ 

19 

READ(1 )NTS 

HPAD 

44 

IF(E0Fd).NE..9)5T0P 

HPAD 

45 

READ(1 )JP 

HPAD 

46 

URITECS, 1 )NTS, JP 

READ C 1 )TST ART. DELT 

HPAD 

47 

HPAD 

48 

15 

READd  )(R(I),I*1,JP> 

HPAD 

49 

RETURN 

HPAD 

59 

1  FORMAT ( 1H9, 1X.24A2, 119 ) 

HPAD 

51 

END 

HPAD 

58 

SUBROUTINE  DAT AS ( KM )  HP AO 

DIMENSION  NTSDC24)  DATAS 

COMMON  R<25000 ), T( 23000 ),P( 23000  ),N0S,LP,LABEl<4  ),NV(6  ),V(6 ),  IQ(S )  HPAD4 

1  ,SBL,  I,TZ, JP, IN, CELT, ID2,H0P< 10 ),ILC,B,SA,SB,DTINE,NJ< 10 ),TSTART,  HPAD 

2  SC,  SKP,  XAN,  XFAC,  YFAC,  TITOO),  ISU,  ICP,  VB,  XI,  NTS(24),ITZ,  CON 

3  JPI(50),IZ, ITN,NOSR,N, IPN0(5 ),SS,KSTS1,KSTS2, VAX, J$U,KU 
DIMENSION  PI(1),PI2<1) 

EOUIUALENCE  <R< 1 ) ,PI C 1 )), (R< 12501 ),PI2( 1 ) ) 

JT-0 
ICN-0 
NSU-1 
REWIND  4 

IF( JSU.EQ.l )G0  TO  10 
N0S1 -NOS 
60  TO  20 
10  NOS1-NOS-1 
20  DO  60  IJ-1.N0S1 

IF(HJ<IJ). £0.0)60  TO  SO 
IF(NSU.CT . 1 )G0  TO  30 
NSU-2 

REA0(4)NT$ 

READ(4)TT,DTINE 
GO  TO  40 
30  READ(4)NTSD 
READ (4 )TTD,DTD 
40  JIP-JTTNJ(IJ) 

JT-JT-H 
READ (4  )N 

READ(4)  (P(IU).IU-JT, JIP) 

JT-JIP 

50  KN-ICN*NJ(IJ ) 

60  CONTINUE 
URITE(3)NTS 
URITEO)KN 
URITE(3)TT,DTINE 
URITE(3)(P(IU),IU-1,KN) 

RET 
END 


1 


c 

c 


SUBROUTINE  PRESS 

CONNON  R(2S000 ),T(2300O),P(23OOO  )»NOS,LP,LABEL(4),NY<6  ), V(6 ),IQ(S ) 

1  ,SBL, I,TZ, JP,  IX,DELT,  Id£,N0P(  10  ),LLCcB,SA,SB,DTIN£,NJ<  10  ),TSTART, 

2  sc,  skp,  xax,  XFAC,  VaC,  TITOS).  isu,  KP.  yb,  XB.  NTS(24),ITZ, 

3  JPI(50),iz,ITN,N0SR,N, IPN0(5  ),SS,KSTS1,KSTS2,VAX,  JSU.KU 
DINENSION  PI(1 ),PI2(1 ) 

EQUIVALENCE  (R( l >,PI( 1 > ), (R( 12501 ),PI2< l ) ) 

CALCULATING  PRESSURE 
N-l 

IZ-IZ+1 

IF  (IZ.GE.ITN)  GO  TO  8* 

30  P(N)-R(I) 

CF-SA+SBX(P(N)-SBL )+SC*(P(N)-S8l)**2 
PINJ-BtCF 

T( N ) -TZ+FLOAT ( IZ >*DTINE 

IF  (IZ.GE.ITN)  GO  TO  40 

IF  (I+IDZ.GT.JP)  GO  TO  40 

N-N+l 

I-I+IDZ 

IZ-IZ+1 

GO  TO  30 

40  IF  (HOP(3).NE.0)  GO  TO  90 
IF  (LP.GT.l >  GO  TO  50 

PRINTING  PRESSURE 

URITE(6,10O)DTXNE 
50  DO  55  NL-1,N,15 
RtN-Nl+14 
IF(HLN.GT.N)NLN-N 

URITE(6,110)T(NL ), (P(NK  ),NK-NL,NLN) 

55  CONTINUE 
CALL  DATAF 

ERROR  PRINT 

IF  (I.GT. JP.AND.LP.EQ.NOS)  GO  TO  50 
GO  TO  90 

60  URITE  (6,120)  I,JP,N 
GO  TO  90 
80  JSU-1 
N-0 

90  RETURN 

100  FORNAT ( 1H0, 4SH  TINE  AND  DATA  POINTS  SAUED  -  TINE  INTERVAL,  NS-, 

*  F8.6) 

110  FORNAT(1X,F10.4, 15F8.2)  _  _  ^ 

120  FORNAT  (1H0,//,2SHI  IS  GREATER  THAN  JP,  I  -, I10,3X,4HJP  -,I10,3HN 
1-,U0) 

END 
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PRESS 

HPAD4 

HPAD 
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HPAD 

HPAD 

HPAD 

PRESS 

PRESS 

PRESS 

HPAD 

CIS 

PRESS 

PRESS 

PRESS 

PRESS 

C17 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

HPAD 

C9 

C10 

C19 

C9 
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C12 

PRESS 

PRESS 

PRESS 

PRESS 

C19 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

PRESS 

HPAD 

HPAD 
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PRESS 
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SUBROUTINE  PLOTS 

COMMON  R ( 2S400 ) » T ( 83000 ) * P < 23000 ) , NOS ,  LP , LABEL (4),NY(6),Y(6),IQ(5) 
1  ,SBL,I,TZ,JP,IX,DELT,ID$,N0P<10>.LLC,B,SA,SB.DTIHE,NJ<10>.TSTART, 

8  sc,  skp7  xax,  xfac.  Vac,  titooi,  isu,  kp,  vb,  xi.  nts<S4>,itz. 

3  JPI (50),  lz,ITN,NOSR,N,  IPN0(5),SS,KSTS1,KSTS2,YAX,  JSU,KU 
DIMENSION  PI ( 1 ),PI3( 1 ) 

EQUIVALENCE  (R( 1 ),PI ( 1 )  >, <R< 12501 ),PI2< 1 > ) 

THIS  SUBROUTINE  PLOTS  REGULAR  PLOTS 

IF  (ISU.GT.l)  GO  TO  10 
ISU-ISW+1 
XPAGE-XAX+3.5 
VPAGE*YAX+3.S 

CALL  PLTBEGIXPACE,  YPAGE, 1 . , 13, LABEL ) 

10  KM-N 

CALL  FIXSCA  (P(  1  ),(CN,YAX,YS,YNI,YMA,DY ) 

CALL  FIXSCA  (TCI ),KH,XAX,XS,XMI,XHA,DX ) 

CALL  PLTSCA<XB,YB,XMI.YNI,XS»YS> 

CALL  PLTAXS<DX»DV,XHI,XMA»YNI,YMA,4> 

TX-XNI-l.SBXS 

CHT*.l 

TY»VNA*1.2iYS 

ENCODE  <2S,140»TIT<1)  )(NTS(IT),IT-1, 18) 

CALL  PlTSVmCHT,TIT(l),0.,TX,TY>  _ 

CALL  LABELA  <DX,DY,XHI,XMA, YNI,YMA,XFAC,YFAC ) 

CALL  PLTUND<XNI,XMA,YMI,YNA> 

CALL  PLTDTS(1,0,T( 1 ),P< 1  ),KN,0) 

CALL  PLTPQE 

IF  CLP.NE.NOS)  GO  TO  130 
IF(N0S.E0.1)GO  TO  130 
REWIND  2 

IF  (JSU.EQ.l)  GO  TO  100 
NOS 1 « NOS 
GO  TO  110 
100  N0S1 *N0S-1 
110  DO  180  IJ-1,N0S1 

IF  (JPI(IJ).EQ.O)  GO  TO  180 
JIP-JPKIJ) 

READ  (8)  (TCIU),IU-1,JIP),(P(IU),IU-1,JIP) 

180  CONTINUE 
130  RETURN 

140  FORMAT ( 12A2, 1H> ) 

END 


PLOTS 

a 

PLOTS 

3 

PLOTS 

4 

HPAD4 

1 

HPAD 

3 

CON 

4 

HPAD 

4 

HPAD 

5 

HPAD 

6 

PL0T8 

7 

PL0T8 

8 

PLOTS 

8 

PL0T8 

18 

PLOTS 

11 

C3 

6 

C3 

7 

C8 

16 

HPAD 

58 

PLOTS 

37 

PLOTS 

38 

CS 

17 

CS 

18 

C4 

3 

CS 

1 

PLOTS 

43 

PLOTS 

45 

CS 

19 

PLOTS 

58 

CS 

28 

C3 

8 

CS 

22 

PLOTS 

68 

HPAD 

59 

PLOTS 

67 

PL0T2 

68 

PLOTS 

69 

PLOTS 

78 

PLOTS 

71 

PLOTS 

72 

PLOTS 

73 

PL0T2 

74 

PLOTS 

75 

PLOTS 

76 

PLOTS 

77 

PLOTS 

78 

HPAD 

68 

PLOTS 

88 

oo 

oo 


1 

5 


I# 

15 


SUBROUTINE  SKIP 

C 

COMMON  R(25###),T(23###),P(23###),N0$,LP,LABEL<4  )«NY(6  ),Y(6),IQ(5  ) 
t  ,  SBL,  I,TZ,  JP,  IX,,DELT#  IDZ,M0P(1#>,LLC,B,SA,S1,DTINE,NJ<  1#  KTSTART, 

2  SC,  SKP,  XAX,  XFAC,  VFAC,  TITO#),  ISU,  ICP,  YB,  XB,  NTS<24),ITZ, 

3  JPI<S#),xZ#  ITM,NOSR,N,  IPN0C5  ),SS,KSTS1 ,KSTS2, YAX, JSU,KU 
DIMENSION  PI( 1 )#PI2( 1 ) 

EQUIVALENCE  <R<  1  ),PI<  1)  >,  (R(  125#1  >,PI2<  1 ) ) 

WRITE (6,1 )IQ 
1##  CALL  DATA  IN 
DO  15#  1-1,5 

IF ( NTS ( I  )«NE.  IQ(  I )  )G0  TO  1## 

15#  CONTINUE 
2##  URITE(6,2 ) 

RETURN 

1  FORMAT ( '  SKIP  TO  CHANNEL  ',3A2> 

2  FORMAT ( *  END  OF  CHANNELS  TO  BE  SKIPPED  * ) 

END 


HP  AD 

61 

HP  AD 

#8 

HP  AD 

#3 

HPAD4 

1 

HP  AD 

3 

COM 

A 

HPAD 

A 

HPAD 

5 

HPAD 

6 

HPAD 

65 

HPAD 

66 

HPAD4 

5 

HPAD4 

6 

HPAD4 

7 

HPAD 

7# 

HPAD 

71 

HPAD 

72 

HPAD 

73 

HPAD 

7  A 
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32S*M33m22288:J*a22832**I«8*S«2 


$IQ(2)  9IQO) 

.339996-94  .2357 IE- 


.91 

-1.99 

-1.99 

-1.99 

.92 

-1.93 

-1.95 

-1.94 

.93 

-1.94 

-1.94 

-1.93 

.99 

-1.92 

-1.92 

-1.99 

.91 

-1.93 

-1.93 

-1.99 

.94 

-1.96 

-1.94 

-.99 

.92 

-1.94 

-1.93 

-1.91 

.94 

-1.95 

-1.94 

-1.91 

.91 

-1.94 

-1.95 

-1.94 

.91 

-1.92 

-1.91 

-1.91 

.91 

-1.92 

-1.91 

-1.99 

.92 

-1.93 

-1.93 

-1.92 

.92 

-1.91 

-1.92 

-1.94 

.99 

-1.94 

-1.94 

-1.99 

.99 

-1.96 

-1.94 

-1.93 

.94 

-1.93 

-1.94 

-1.97 

.95 

-1.94 

-1.95 

-1.95 

.71 

-.69 

-.69 

-.99 

.72 

-.73 

-.72 

-.71 

.59 

-.71 

-.79 

-.69 

.71 

-.7a 

-.72 

-.71 

.71 

-.71 

-.72 

-.71 

.71 

-.69 

-.69 

-.71 

.69 

-.69 

-.69 

-.99 

.69 

-.69 

-.79 

-.7* 

.69 

-.79 

-.72 

-.74 

.76 

-.76 

-.74 

-.73 

.74 

-.73 

-.72 

-.78 

.75 

-.75 

-.76 

-.79 

.97 

-.99 

-.96 

-.93 

.99 

-1.92 

-1.94 

-1.93 

.97 

-1.99 

-1.99 

-1.95 

.99 

-1.11 

-1.11 

-1.99 

.96 

-1.99 

-1.95 

-1.91 

.99 

-.92 

-.93 

-.92 

.96 

-.97 

-.99 

-.29 

.97 

-.97 

-.99 

-.22 

.95 

-.99 

-.92 

-.91 

.99 

-.91 

-.93 

-.22 

.99 

-.92 

-.94 

-.24 

.79 

-.77 

-.76 

-.79 
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TABLE  El.  INPUT  CARD  IMAGE 


Variable 

Format 

Column 

Description 

IPNO( 2) 

ii 

4 

Sets  IDZ  and  IX  as  specified  in 
Table  E2. 

IDZ 

The  increment  for  choosing  the 

input  data  to  be  saved 

(e.g.,  5=save  every  fifth  sample.) 

IX 

The  number  of  samples  to  skip 
after  the  fiducial  mark  to  reach 
the  data  which  are  to  be  processed 

IPNO( 3) 

ii 

5 

Sets  KP,  M0P(3)  and  M0P(4)  as 
specified  in  Table  E2 • 

KP 

Plotting  option: 

0  =  no  plot, 

1  =  special  plot  to  be  programmed 

by  user, 

2  «  normal  page-size  plot. 

MOP( 3) ,M0P(4) 

Choices  for  integrating:  none, 
once,  or  twice* 

B 

F10.3 

6-15 

Calibration  constant  for  the  gage 
in  engineering  units  per 
calibration  step* 

IXS 

110 

16-25 

0  =  no  change  to  the  value 
set  in  IPN0( 2) ; 

0  *  change  the  value  set  in 

IPN0( 2)  to  this  value. 

DLTM 

F10.0 

26-35 

Time  adjustment  to  be  subtracted 
from  the  value  of  TSTART. 

IDZS 

15 

36-40 

0  -  no  change  to  the  value  set 
in  IPN0( 2) ; 

0  t  change  the  value  set  in 
IPN0(2)  to  this  value. 

NLIST 

15 

41-45 

0  «  NAMELIST  not  used; 

1  =  NAMELIST  used  (see  Table  E3.) 

IFID 

15 

46-50 

0  =  fiducial  mark  is  present; 

1  *  no  fiducial  mark. 

NLK 

15 

51-55 

Number  of  contiguous  data  sets  to 
be  processed  the  same  way. 
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TABLE  E2.  PARAMETER  SETTINGS 


Parameter 


Value  Description 


IPNO( 2)*  0 


1 

2 

3 

4 
3 
6 

7 

8 

9 


Skip  this  data  set  if 

IPNO( 3)  =  0  also;  otherwise, 

an  error  message  is 

printed  and  processing  halts 

IDZ  =  1  and  IX  -  500 

IDZ  =  1  and  IX  =  2000 

IDZ  =  1  and  IX  =  3000 

IDZ  -  2  and  IX  =  500 

IDZ  =  2  and  IX  =  2000 

IDZ  =  2  and  IX  =  3000 

IDZ  =  5  and  IX  =  500 

IDZ  -  5  and  IX  =  2000 

IDZ  =  5  and  IX  =  3000 


IPN0( 3)** 


0 


1 

2 

3 

4 

5 

6 

7 

8 

9 


Skip  this  data  set  if 
IPN0(2)  =  0  also;  otherwise, 
an  error  message  is  printed 


and  processing  halts 
KP  ■=  0,  MOP ( 3 )  -  0*** 
KP  =  0,  M0P( 3)  =  2, 

KP  =  0,  MOP( 3)  =  2, 

KP  =  1,  MOP( 3)  =  0, 

KP  =  1,  M0P( 3)  =  2, 

KP  =  1,  M0P( 3)  =  2, 

KP  =  2,  M0P( 3)  =  0, 

KP  *  2,  M0P( 3)  =  2, 

KP  =  2,  MOP( 3)  =  2, 


and  MOP( 4)  =  0*** 
M0P(4)  =  0 
M0P(4)  =  2 
MOP(4)  -  0 
M0P(4)  =  0 
M0P(4)  =  2 
M0P(4)  ■  0 
M0P(4)  =  0 
MOP(4)  =  2 


*  If  these  choices  for  IDZ  and  IX  are  unacceptable,  choose  the  option 

closest  to  the  desired  value  and  insert  the  correct  values  as  IDZS  and/or 
IXS  on  the  input  card  image. 

**  If  IPN0( 2)  =  0  and  IPN0(3)  =  0,  the  array  IQ  is  read  in  from  the  next  card 
image  with  format  5A2.  The  array  IQ  must  exactly  match  NAMEF,  as  defined 
in  the  data  file  created  in  ADTAPE2 ,  of  the  next  data  set  to  be 
processed.  If  IQ  and  NAMEF  do  not  match,  the  search  will  continue  to  the 
end-of-file  marker. 

***  0  *  no ;  2  =  yes.  For  one  integration  M0P(3)  =  2  and  M0P(4)  =  0.  For  two 
integrations  M0P(3)  =  2  and  M0P(4)  =  2. 
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TABLE  E3.  NAMELIST  VARIABLES 


Variable 

Default  Value 

Description 

IBSE 

500 

Number  of  samples  to  skip  after  the  fiducial 
mark  before  starting  to  sample  the  baseline. 

ITM 

23000 

Maximum  number  of  data  samples  to  be 
converted  to  engineering  units. 

ITZ 

950 

Number  of  samples  to  skip  after  sampling 
the  last  calibration  step  before  starting 
to  search  for  the  fiducial  mark. 

MOP 

Control  variable: 

0  -  end  of  processing, 

1  ■  conversion  to  engineering  units, 

2  -  integration, 

3  ■  plotting. 

NS 

6 

Number  of  calibration  steps  to  be  used 
in  the  least  squares  regression. 

NX 

40 

Number  of  samples  to  be  averaged  for  each 

calibration  step. 


NY 

16,100,240, 

360,480,600 

Sample  indices  marking  the  position  on 
each  calibration  step  to  start 
processing  data. 

SKP 

150. 

Criterion  for  eliminating  spikes  in 
calibration  steps. 

SS 

1. 

Step  size  of  each  calibration  step, 
usually,  in  ohms  or  volts. 

XAX 

7.5 

Abscissa  graph  size  in  inches. 

XFAC 

1. 

Scale  factor  for  the  numeric  label  on  the 
x-axis  (e.g.,  if  XFAC  *  .001  and  the 
number  to  be  printed  is  10000,  a  10 
will  be  printed.) 

YAX 

5. 

Ordinate  graph  size  in  inches. 

YFAC 

1. 

Scale  factor  for  the  numeric  label  on  the 

y-axls,  similar  to  XFAC. 
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DISTRIBUTION  LIST 


No.  of 
Copies 


Organization 


No.  of 
Copies 


Organization 


12 


1 


1 
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USA  ARRADCOM 
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Commander 
USA  ARRADCOM 
ATTN:  DRDAR-SCM 


DRDAR-SCM, 

Dr. 

E. 

Bloore 

DRDAR-SCM, 

Mr. 

J. 

Mulherin 

DRDAR-SCS , 

Mr. 

B. 

Brodman 

DRDAR-SCS , 

Dr. 

T. 

Hung 

DRDAR-SCA , 

Mr. 

S. 

Jacobson 

DRDAR-SCA , 

Mr. 

w. 

Gadomski 

DRDAR-SCA, 

Mr. 

E. 

Malatesta 

Dover,  NJ  07801 


1  Commander  3  Commander 

US  Army  Materiel  Development  USA  ARRADCOM 
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l 


2 


6 


l 


Commander  1+ 

USA  ARRADCOM 

ATTN:  DRDAR-LC ,  J.T.  Frasier 

Dover,  NJ  07001 

Commander 
USA  ARRADCOM 
ATTN :  DRDAR-TSS 

Dover,  NJ  07801  1 

Commander 
USA  ARRADCOM 

ATTN:  DRDAR-LCU ,  Mr.  E.  Barrieres 

DRDAR-LCU ,  Mr.  R.  Davitt 
DRDAR-LCU-M ,  Mr.  D.  Robertson  2 
DRDAR-LCU-M,  Mr.  J.  Sikra 
DRDAR-LCU-M,  Mr.  M.  Weinstock 
DRDAR-LCA,  Mr.  C.  Larson 
Dover ,  NJ  07801 

Commander 
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